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We use variational perturbation theory to calculate various universal amplitude ratios above and 

below Tc in minimally subtracted (/i^-theory with A*" components in three dimensions. In order to 

c 2 ^ ' best exhibit the method as a powerful alternative to Borel resummation techniques, we consider 

o ■ 

' only to two- and three-loops expressions where our results are analytic expressions. For the critical 

^ ■ exponents, we also extend existing analytic expressions for two loops to three loops. 

o 
o 

(N 



I. INTRODUCTION 



' Recently, quantum mechanical variational perturbation theory W has been successfully extended to quantum field 
theory, where it has proven to be a powerful tool for determining critical exponents in three iQ-Q as well as in 4 — e 



■ dimensions [g|J^ . The purpose of this paper it to apply this theory to amplitude ratios which can be measured exper- 

o : . 

: imentally. Their perturbation expansions suffer from the same asymptotic nature as those of the critical exponents, 

^ . thus requiring delicate resummation procedures. These have been the subject of numerous studies, of which we can 

^ only mention a few, by various groups. There are two main approaches followed by various authors which we shall 

Q ' divide according to their method into a Paris school and a Parisi school. 
O ■ 



^ ' The Paris school follows Wilsons ideas |j7|j8| by considering epsilon expansions in D — 4 — e dimensions, making 
•I— I , 

use of the fact that in the upper critical dimension D^p = 4 the theory is scale invariant. The results are at first 

H ■ 

power series in the renormalizcd coupling constant g. For small e, the coupling constant goes, in the critical limit of 
vanishing mass, to a stable infrared fixed point g — » 5*, where scaling laws are found The position of the fixed 
point is found as a power series in e which makes critical exponents and amplitude functions likewise power series in 
e. These series diverge. The large-order behavior ^0 11| suggests that these series are Borel summablc [^^. The 



exact e-expansions of the critical exponents arc known up to five loops ||lJ,|T5|. They have been resummcd with the 
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help of Borel transformations and analytic mapping methods in Refs. ]16|-|18| 



The Parisi school follows Ref. |19| in studying perturbation expansions directly in 13 = 3 dimensions [gO-25|. In 
the original works, renormalization conditions are used according to which renormalized correlation functions should 
behave for small momenta like G{p) ~ {p^ + m^)^^. Recently, these normalization conditions have been replaced 
by dimensional regularization near _D = 3 to remove divergences (see for instance p3| which uses a regularization in 
= 3 — 8 dimensions). Contrary to the e-expansions around D^p = 4, the system is no treated near the dimension 
of naive scale-invariance, and the scaling properties are no longer obvious order by order in g. In addition, singular 
terms violating Griffith's analyticity are introduced which show up by amplitudes having unpleasant logarithmic 
dependences on the coupling constant. 



The universal amplitude ratios were first discussed in 1 26 1 in the context of Wilson's renormalization group approach, 



and by Bervillier p7[ within the field theoretic approach developed in The experimentally most easily accessible 
amplitude ratios are formed from the amplitudes of the leading power behaviors of various physical quantities in 
T — Tc above and below the critical temperature Tc- A typical example, and one of the best measured amplitude 
ratios, is for the specific heat of superfluid helium above and below Tc- It was obtained in a zero-gravity experiment 
by Lipa et al. who parameterized the specific heat as follows (we use the second of the references quoted in p8[): 

C^=A^\t\-'^[l + D\t\^ + E\t\^^) + B, t = T/T,-l, (1) 

with a = -0.01056 ± 0.0004, A = 0.5, A+/A- = 1.0442 ± 0.001, A" = 525.03, £> = -0.00687,^; = 0.2152 and 
B = 538.55 (J/mol K). This parametrization is an approximation to the Wegner expansion form 

F = F±\t\>^{l + aoMt\^° + ao,2\t\'^" + aoMA'^" + ■■■ 

+ ai,i|i|^i+ai,2|ip'^^+ai,3|tP'^^+---) (2) 

with X ^ combination of critical exponents and F± denoting the leading amplitude above and below Tc, respectively. 
Compared to this general Wegner expansion, higher powers in Aq = A have been neglected in (^, as well as daughter 
powers A.i,i > 1. This will be also the case in the present work, where we shall take into account only one exponent 
A, related to uj by the relation A = uny. 

Further amplitude ratios are formed from the amplitudes a^j of the nonleading power behaviors in T — Tc, the 
so-called confluent terms or crossover functions, these also being universal quantities p^ , pO[ . They arc known up to 
three loops. None of them will be examined here. 



Apart from critical exponents and amplitude ratios, experimental observations show that the equation of state and 
the free energy have a simple scaling form of the Widom type, whose field-theoretic explanation can be found in 
various textbooks For example, the free energy of a system with magnetization Mb may be represented 



the critical temperature. The scaling equation of state has been calculated in e-expansions to order for general 
0(-/V)-symmetry Q and to order e' for the Ising model {N — 1) 



Amplitude ratios relate the properties of the disordered phase, which are easy to calculate, to those of the ordered 
phase, which arc much harder to derive. Several methods have been proposed to connect the two phases. One of them 
is due to Bagnuls and Bervillier [ p4[ , and was applied further in p5[ | . A similar procedure was followed in p^ , |23t for 
the amplitude ratio of correlation lengths, which had been omitted by Bagnuls and Bervillier. Calculations in three 



three-loop results). Such analytic studies are important since they offer insight into the nonanalyticity with respect 
to the coupling constant. The amplitude ratio found in is restricted to the Ising case = 1. The same is true 

for 1^^, which includes all diagrams up to three loops. 

All power series are divergent and require resummation. Numerically, this has been done for the Ising model in 
pl| to five loops for the critical exponents, various amplitude ratios, and the equation of state. Reference |^l| also 
contains comparisons between the results of different groups (both for D ~ i and D = 4 — e), with experiments 
and with high-temperature series. For the most up-to-date work, see [|l8| , which besides the critical exponents and 
amplitude ratios for the Ising model gives also the critical exponents for general O(A^) symmetry. 

An other approach has been followed by Dohm and collaborators in Aachen who proposed to use an 

analytic renormalization scheme in the form of minimal subtraction when working in D ~ i dimensions. The use of 
the minimal subtraction scheme in field theories at fixed dimensions 2 < _D < 4 has one important advantage: the 
renormalization constants are the same in both the symmetric phase with T > Tc and the ordered phase with T < T^. 
The renormalization constants are power series in the renormalized coupling constant with coefficients which are poles 
in e up to the given order of the perturbative series 




A. Perturbative calculation of amplitude ratios 




||,|| for analytic 



L 





(3) 



The most important property of this scheme is that the mass does not enter exphcitly the expansions, which can 
therefore be used on both side of Tc. Since the critical exponents are related to the renormalization constants, the 
mass independence of the Zi implies a clear decomposition of the correlation functions into amplitude functions and 
power parts. Working in three dimensions, there is a prize to pay: logarithmic singularities in the coupling constant. 
They can be removed using suitable length scales. This may be the physical length scale ^+ above Tc, and an other 
length scale ^_ related, in the critical regime, to the longitudinal mass below T^. Since they are not exactly equals, 
the Aachen group call the length scale a pseudolength. A precise definition of ^_ has been given in ||3^. With 
different collaborators, Dohm has applied this scheme to derive various critical exponents and renormalization-group 
functions above Tc |35| , to calculate the heat capacity, the order parameter and the superfluid density (both above 



and below Tc)^ as well as some useful universal combination of observable quantities |36|. So far, these works have 
been limited to low orders. The Ising model is the simplest system, since it contains no massless Goldstone modes 
which cause extra infrared singularities at intermediate stages of perturbative calculations of the thermodynamical 
quantities on the coexistence curve where the external magnetic field vanishes. The infrared singularities are the 

Is^JstII to two 



reason why the analytical equation of state and amplitude functions below Tc have been restricted 1 27 
loops for general N . These extra infrared singularities, which cancel at the end of the calculations, are caused by 
the physical singularities of the transverse susceptibility. Being physical, they remain at the end. Due to these 
difficulties, numerical studies up to five loops below Tc, with accurate Borel resummation are available only for the 
Ising case pl| , p5|p^j39| . Only analytic three- loop calculations for the thermodynamic quantities below Tc have become 
recently available for the general 0(A^)-system Based on these, calculations in which some contributions were 
evaluated up to five loops were done for amplitude ratios at iV = 2 and iV = 3 |^^, proceeding as follow: Amplitude 
functions for the heat capacity were calculated using the three-loop result of ji^ and five-loop results for the vacuum 
renormalization constant [^9|j4l| ] and the critical exponent a. For a, use has been made of the values given in |l^] 
for N — 1, of the value given in the first of Ref. |28j for N ^ 2 (this being the initial result of the space shuttle 
experiment, which was subsequently corrected), and of the value given in Ref. for iV = 3. Since then, the works 
LSljf^] have appeared and seem to be the best available references concerning resummcd data. Although this is 
not the main subject of this paper, it is interesting to see in which way the new values of a affect the amplitude ratios 



of the heat capacity given in ||39|]. This will be done in Section [V. 

In the following, we shall calculate amplitude ratios with the help of Kleinert's variational perturbation theory 
L3| . To exhibit the method most clearly, we shall base our study on analytical results only. This will restrict us 



to the level of three loops. Working at such low orders, the accuracy of our resummed values cannot compete with 
some existing five-loop calculations. For this reason, we shall not include nor discuss error bars in the final results. 

To illustrate the method of variational perturbation theory, we shall first show how to obtain analytic expressions 
for the critical exponents, thus extending an earlier two-loop analytic calculation in Ref. After this, we apply the 
procedure to amplitude ratios of various experimental quantities. The critical exponents are computed directly from 
the renormalization constants of the theory. In the minimal subtraction scheme, the renormalization constants have 
only pole terms in e. For the amplitude functions, this is no longer true: in a D = 4 — e approach, they have to be 
expanded in e. For this reason it is not a priori clear at which level the variational method has to be applied. For 
the purpose of showing the power of the method to resum amplitude ratios, it is then better to calculate amplitude 
ratios in three dimensions. A resummation of amplitude ratios within the e-expansion method is postpone to a later 
publication We shall also consider only the expansions of the Aachen group, especially their analytical two- loop 
and three-loop expansions. As a bonus, since the renormalization constants are the same (apart for trivial 
coefficients coming from the respective conventions) in the minimal subtraction scheme in _D = 4 — e dimensions and 
in fixed D = S dimensions, the critical exponents will be the same in variational perturbation theory. This will be 
shown explicitly below. 

The paper is organized as follow. In Section ||, we define the model and the conventions. In Section III, we briefly 



review the strong-coupling approach and apply it to the evaluation of the critical exponents at the level of two and 
three loops, extending the results of Ref. ||]. Section IV is the main part of this paper, where we show how the 



strong-coupling limit of various amplitudes and amplitude ratios are determined. In Section IV F , we use the latest 
available value for the exponents a and v P^ Jl^ , |l8|] to calculate the amplitude ratio of the heat capacity and the 
univeral combination Rc (constructed from the leading amplitudes of the heat capacity, the order parameter and the 
susceptibility above Tc), ior N = 0, ■ ■ ■ ,A, and to calculate the amplitude ratio of the susceptibilities in the Ising model 
(N = 1). Finally, we draw our conclusion in Section For completeness, we have added an Appendix containing 
all formulas taken from other publications, and calculations related to them. 



II. MODEL AND CONVENTIONS 

The critical behavior of many different physical systems can be described by an 0(A^)-symmetric (/)^-theory. In 
particular, the case = describes polymers, = 1 the Ising transition (a universality class which comprises 
binary fluids, liquid- vapor transitions and antiferromagnets) , N — 2 the superfluid Helium transition, A^ = 3 isotropic 



magnets (transition of the Heisenberg type), and = 4 phase transition of Higgs fields at finite temperature. In the 
presence of an external field hs , the field energy is given by the Ginzburg-Landau functional 



(4) 



To facilitate comparisons with the results of the Aachen group we use the same normalizations. The fields 

(j)B and the external magnetic field Hb have TV components, is the bare coupling constant, and ro a bare mass 
term, to be specified later. The integrals are evaluated in dimensional regularization. In dimension _D = 3, 0^-theory 
is superrenormalizable. This means that only a finite number of counterterms have to be added in order to make 
observables finite. More economically, the divergencies can be removed by a shift of the mass term and reexpanding in 
''0 ~ foc, where tqc is the critical value of tq. In e-expansions, tqc vanishes. Near the critical temperature, rg behaves 
like roc + aot, where t is the reduced temperature (T — Tc)/Tc. When working near D = 3 dimensions, it is possible 
to use a simplified shift 6rQ that only contains the D — 3 pole of tqc (and not the poles at Z?; > 3 with / = 3, 4, 5, • • •, 
where Di = 4 — 2/1). For convenience, wc write the differences as a new mass term: ro — roc = ni^ and ro — (5ro — m'^B- 
In this way, we arrive to a new bare theory, with a mass term rre'^ which may be considered as the physical square 
mass of the theory. The introduction of the mass m'g makes the theory finite. It has however to be distinguished from 
the mass, field and coupling constant renormalization which still has to be performed: this latter renormalization, 
related to the introduction of the renormalization constants Zi, is nothing else than a change of variables reflecting 
the fundamental scale-invariance hypothesis of the renormalization group approach. The distinction between the two 
steps - making the theory finite and renormalizing - is irrelevant in D = 4 — e dimensions because roc = at e = 0: 
Finiteness of the theory and the renormalization program are more intimately related than in D = 3 dimensions. For 
a thorough discussion of the difference between the renormalization in D = 4 — e and fixed D ~ 3 dimensions, see 
[||,|5|, in particular p. 7215 in Q. 

Within the minimal renormalization scheme, the renormalization constants Z, which are introduced to remove the 
poles at -D = 4, are given by 

mB=m -—, (5) 
Adub = iJ-^-^^u, (6) 

= (7) 



the quantities on the right-hand-side being the renormalized ones. In Eq. (g), /i is an arbitrary reference mass scale 
and 



Ad = r(l + e/2)r(l - e/2)S'r,, with So 



2^^/2 



7 

(8) 



V{D/2){2-k)D 

is a convenient geometric factor. The number is equal to SdH'^t^)^ where Sd is the surface of a sphere in D- 
dimensions. Since An goes to Sd when Z) — > 4, the renormahzation constants have the same form in Z) = 3 as 
in Z) = 4 — e, and the resummation for the critical exponents is identical for the two approaches. This will be made 
clear below. For the amplitude calculations, however, things are different: If the expansions are truncated at some 
order, they turn out to depend on the difference between Ad and Sd- Rather than saying that the normalization of 



Ad is a matter of convenience to simplify the D-dependence of lower order results 35 1, we shall see that the use of 
the geometric factor improves low-order results: For example, the one-loop expansion of the amplitude function 
for the order parameter is identical to the zero- loop order js^. 

With these conventions and notations, the rcnormalization constants in minimal subtraction are given up to three 
loops by m-lll] 



Z™2 — 1 



4(7V-f-2) 



+ 8(iV-K2) 



2(7V + 5) 



4-8(iV + 2) 



Zu = l + — — + 16 



d,{N + h){N + &) 4(lliV-h50) 31iV-H230 



2 



5A^ + 22 



24(iV + 8)3 16(A^ + 8)(17iV + 76) ^ 96C(3)(57V + 22) + ibN"^ + 942A^ + 2992' 
4(iV + 2) o 8_, _ /4 1 



-{N + 2){N + 8) \ 



(9) 

(10) 
(11) 



They are related to that in Ref. [|3j by the replacement u ^ 5/12. This factor comes from the different coefficient of 
the coupling term u g/Al in (Q) and the fact that a factor l/(47r)^ is absorbed in the definition of g in ||l^, whereas 
a factor Ad=a — l/(87r^) is included here. 

These rcnormalization constants serve to calculate all critical exponents including the exponent lo which charac- 



terizes the approach to scaling. This is the subject of Section [I] in which we illustrate the working of variational 
perturbation theory. 



III. EXACT CRITICAL EXPONENTS UP TO THREE LOOPS 

Variational perturbation theory has been developed for the calculation of critical exponents in |^ and |^ in D = 3 
and D = A — e dimensions, respectively. A review can be found in the textbook JlSf . So we need to recall here only 
the main steps of the procedure. 



Let fhiuB) be the partial sum of order _L of a power series 

L 

f^fL{uB)^Y.f^^B- (12) 

1=0 

In the present context, 

ub = UBfJ-^^Ao (13) 
with D = 3 and e = 1, i.e., ub — ws/(47r/i). The mass scale ^ will be specified later. As seen from Eq. (0), this scale 



leads to a dimensionless coupling constant ub- We assume that in Eq. (12), the ultraviolet divergencies have been 
removed. In _D = 3 dimensions, this is achieved by working with rrig instead of tq. However, tqc is a nonperturbative 
quantity in three dimensions, and working with or m'^ generates nonanalyticities due to the presence of logarithms 
of the coupling constant. These will be removed by the introduction of the correlation length C+ above Tc and of the 
length below Tc, see [Q. The mass scale fj, will be identified with the inverse of these correlation lengths C±^ in 
the two phases. Since the correlations lengths go to infinity like \t\^'^ as the critical point is approached, the series 
have to be evaluated in the limit of an infinite dimensionless bare coupling constant ub- In the rcnormalization group 
approach, this regime is studied by mapping the expressions into a regime of finite renormalized quantities using the 
renormalization constants (^-(0). If we can find directly the strong-coupling limit, this rcnormalization is avoidable. 
To understand this, consider the relation between the renormalized and the bare coupling constant at the one-loop 
order u ~ ubIJ^ '"' — c/ e{uBlJ''^Y , where c is a constant. At the critical point, /i ^ 0, or ub oo, and the series 
expansion breaks down. If we sum a ladder of loop diagrams, we obtain l/u = l/(uB/i^'^) + c/e. Now critical theory 
can easily be reached to give a renormalized u* — e/c. A strong-coupling expansion in the bare coupling will turn out 
to give the same result. From our point of view, the renormalization group approach is simply a specific procedure of 
evaluating power series in the strong-coupling limit. 

In _D = 4 — e dimensions, the situation is slightly more involved since renormalization is also necessary to obtain 
UV- finite quantities, the mass shift tq — rpc not being sufficient for this goal as in the superrenormalizable case D — 3, 
since rpc = as e 0. As far as this paper is concerned, we shall make use of the fact that D ~ 3 and D = 4 — e 
dimensions series expansions in term of renormalized quantities are available in the literature. These will be converted 
back to bare expansion, using the inverse of Eqs. For D — 3 dimensions, this expresses all physical quantities 

in powers of ub/ij-- The mass scale /i is identified with in the disordered or ordered phase, respectively. In 
_D = 4 — e dimensions, the critical theory is obtained by identifying —t m with the renormalized mass m in the 
disordered phase. In a subsequent publication [Eq], we will show how to perform directly a calculation in term of 



UV-finite bare quantities in _D = 4 — e. In this way, the renormahzation procedure is superfluous, our sole problem 
being the evaluation of the expansions in the limit of infinite coupling constant. 
Inverting Eq. we have the expansion 

"2(iV + 8)2 3(3iV+14)" 



_ I 4(iV + 8) _ 
u — UbS 1 ■f^-B 



+ 



- 8 



8{N + 8f 32(7V + 8)(3iV+ 14) 96C(3)(57V + 22) + 33iV2 + 922A^ + 2960' 



3e 



(14) 



The expansion ( p^ ) has the same strong-coupling limit in _D = 3 and D ~ 4 — e dimensions, and it does not matter 
that /i = C±^ for D = 3 or ji = m for D = i — e since both quantities go to zero in the critical limit with the same 
power With relation (|l|) between u and ub , we obtain the bare coupling expansion of the renormalized square 

mass and fields; 



2 _ ry-l 2 2 

in = nip, = rrir. 



1- ^^^^^^^ MB+4(7V + 2) 



4(7V + 5) 5 

^2 +- 



l&{N + 2) 



4(7V + 5)(Ar + 6) ^ 53iV + 274 _^ (57V + 37) 



3e2 



■ 2(iV + 2)_2 

1 + Ug 



(15) 
(16) 



These two expressions are sufficient to calculate the critical exponents v and 7 and, via scaling relations, all other 
exponents. Note that the value of the renormalized coupling constant at the critical point u* is not needed to obtain 



V and 7. The expansion (14) is however useful for obtaining an accurate exponent lo of the approach to scaling. It was 
pointed out in ]4^ that lo can also be deduced from the expansions of v and 7. However, to reach the same accuracy, 
this requires always one more loop compared to the loop order we are interested in. For this reason, we shall take the 
advantage of Eq. (M), whose three-loop order contains all necessary information to get uj to that given order. 



A. Method 



Starting from Eq. (p^), we follow jU^,]!^ to write its strong-coupling limit as 



/l(ub ^ 00) = optji^ 



L L-i 
1=0 j=0 



—iuj/e 
j 



{-ly 



(17) 



The symbol opt^^ denotes optimization with respect of ub- This expression holds provided it yields a nonzero 
constant. This limit will be denoted by /*: 



f(uB ^ 00) = /* + CoUg"^^' + Oiug^"^'^), 



(18) 
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where cq is a constant. The optimalization is supposed to make / depend minimally on iiB- In practice, this amounts 
to taking the first derivative to zero (odd orders) or, when it yields complex results, to taking the second derivative 
to zero and selecting turning points. 

After having determined the optimum at various order L, it is still necessary to extrapolate the result to infinite 
order i — > oo. The general large-L behavior of the strong-coupling limit has been derived from an analysis in the 
complex plane in : 

/£«r + ciexp(-c2Li-"), (19) 
with constants ci and C2 > 0. Knowing this behavior, a graphical extrapolation procedure may be used to find 

r.o = r. 

To apply the above algorithm to critical exponents, we proceed as follows: Let Wl be a function obtained from 
perturbation theory. It has an expansion 

L 

Wl{ub) ^Y.^^'^B- (20) 

Suppose that we also know this function has a leading power behavior uF^'' for large ub- The power p/q is given by 
a logarithmic derivative 

q d log Ub 

The right-hand-side is a power series representation of a function of the type (p^, with p/q being /* and the approach 
to /* in the form of powers u^^'^ . Equation (^l|) will be used later for the determination of the critical exponents. 



If the series (^0|) goes to a constant in the strong-coupling limit, the exponent p is vanishing, and we are left with 

1^ = 0. (22, 

a log UB 

This equation can be solved for q, i.e., for uj. Note that (|2l| ) strictly holds for p ^ 0. However, it can be shown that 
this equation may be used also for p — 0, i.e., that ( p2| ) is a consistent equation for functions which go to a constant 
in the strong-coupling limit. This is explained in Appendix ^ In the following, we shall directly use (^) and ( |2^ ) 
for two- and three-loop expansions where everything can be calculated analytically. We give below the associated 
formulas resulting from Eq. ([l7|). Setting p = 1 + e/u, we find for L = 2: 

11=2 = opt^^ (/o + hpUB + f2U%) = /o - (23) 

4 J2 
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while the three-loop results L = 3 leads to 

/£=3 = 0Pt^« (/o + fiUB + hul + hu%) = /" - (l - l^) + (1 - (24) 

where /i = /ip(p + l)/2, /2 = /2(2p — 1), r = ■\/l — 3/1/3//I. If the square root is imaginary, the optimal value is 
given by the unique turning point. Practically, and this is a virtue of the analytic result, this square root is always 
imaginary for Z? = 3, at least as for the exponent ld. The turning point condition leads to 

3/3 h 

i.e., to same expression as Eq. (^J), but with r ~ Q. In the case D = A — e with e ^ 0, r is real. However, for uj the 
e-expansion of r produces higher orders in e than the three-loop approximation admits. Then, in both _D = 3 and the 
e-expansion, ( |25| ) is the relevant equation. A word of caution is nevertheless necessary: The positive root r of 

u*^ = I^[-l±r) (26) 

has to be chosen in order to match the three-loop result with the two-loop one in the limit /a — > 0. Doing so, it must 
be assumed that /2 and /i are nonvanishing. When optimizing with /2 = 0, it is immediate to show that if /1/3 > 0, 
then the optimum corresponds to u|j(/2 — > 0) = and /£^3 = /o- The other possibility, /i = 0, is also interesting 
since it occurs in the determination of the exponent 77. It can be verified that /i = implies taking the negative root 
r — —1, so that ^ 0) = -2/2/(3/3) and /s = /o + 4/|/(27/|). This possibility has not been discussed in the 

previous works PJsljl^. 



B. Critical exponents 

After the introduction to the resummation method to be used in this work, we can now turn to the actual deter- 
mination of the critical exponents. We start from the definitions within the conventional renormalization formalism 
of the functions 



7m 



ji dm? 



vn? 


M 


d 


lo 






u 






B 



(27) 
(28) 
(29) 
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which, in the critical regime 0, render the critical exponents rjm = and ry = 7^ if the first two equations 

are calculated at the fixed point u* determined by the zero of the third function Pu- The derivative of at u* is the 
critical exponent of the approach to scaling lu — d(3u/du\^,. 

Using the relation between the bare coupling constant ub and the reduced one ub given in Eq. (p^, Eqs. ( |27| ) 
and ( |2^ ) become 

d , rn? d , 

'n^ = -'^T^ ^log— =="*^Ti ^log^r > (30) 

d log u B nig d log u b 

77 ==£71 ^log^r^^e— ^logZ, ' , (31 

dlogUB (I)b dlogUB ^ 

where the renormalization constants Z^^ and have been explicitly given up to three loops in Eqs. ( p^ ) and ([l^), 

respectively. The associated power series expansion in ub of the exponents rjm and 77 will now be treated with the 
help of the formalism described in the previous section, up to two and three loops. 



C. Critical exponents from two- loop expansions 



In order to calculate the two-loop expansions in the critical strong-coupling limit, we need to know uj to this order. 
This will be calculated from Eq. (p^). Dividing this series by ub, we know that the leading power behaviour as 
Ub ^ cxd is —1 since u is supposed to go to the constant value u*: uug^\^^^^ — u*Ug^. Calculating the logarithmic 
derivative of (H) and expanding up to second order in ub, we have 



d , u 
-J] — ^log — 
d log Ms Ub 



-A{N + 8) 



Ub + 16 



(TV -I- 8)2 3(3A^+14) 



u%. 



We now apply formula (A4). Combining with (E3), we identify 



[-4(iV + 8)/e]' 



4 16 [{N + 8)2/e2 + 3(37v + 14)e] ' 



(32) 



(33) 



I.e., 



from which we can deduce lo: 



4 



3e 



3iV-h 14 
'(A^ + 8)2 



(34) 



(35) 



p-1 -1 + 2^1 + 3e(3iV + 14)/(7V + 8)2 ' 
It is identical to the result obtained in |^ . As a check of (|35| ) , we verify that it reproduces the well-known e-expansion 



uj, = e- 3e2(3iV + 14)/(iV + 8)^ 
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(36) 



We refer the reader to |§|,|T3| for plots of the function ( p5|) as e goes from to 1, and for a comparison with the 



unresummed e-expansion. The strong-coupHng Umit of uj may also be calculated from ( Al) with an analytic expression 
different from (|35|), although numerically they are practically the same, and they certainly have the same e-expansion 



This determination of lj illustrates what we said in Section ||, that in the minimal renormalization scheme the 
critical exponents lead to identical results in Z) = 3 and D = A — e dimensions. This will also be true for the 
critical exponents to be calculated in the sequel [Q. For this reason, we shall always keep trace of e to facilitate the 
comparison, although our work is in Z? = 3 dimensions. Only for amplitude ratios to be calculated later will such a 
comparison be impossible and e be set equal to 1 everywhere. 



Knowing w, we can now determine the exponents r] and rjm- According to (30) and (|3l|), we take the logarithmic 
derivative of (|l5|) and (|l^), reexpand the results up to the second order in -u^, and obtain 

"2(A^ + 8) 



ry™ = 4(iV + 2)uB - 8(A^ + 2) 



'■Bi 



= 8{N + 2)i 



■B- 



(37) 
(38) 



Evaluating rjm in the strong-coupling limit in the same way as w, i.e., following the algorithm (|T^), we find 



Vrn = 



m+2)Y 



{N + 2) 



,3(37V + 14) 
(iV-l-8)2 



(39) 



4 8{N + 2) [2(iV -t- 8)/e + 5] {N + 8) + 5e/2 
For 77, the situation is less clear. In ||^,|^, it was argued that the two-loop result cannot be computed from Eq. ( |3^ ) 
since no linear term in ub is present. A direct application of the resummation algorithm would give an optimum 
u*g = 0, then a value t] — at two-loop order. This does not lead to the correct e-expansion, according to which the 
exponent start with e^, i.e., with a non- vanishing two- loop contribution. To apply variational perturbation theory, it 
is necessary to modify the procedure. In Ref. ||^, this was done by considering a different critical exponent 



7 = 1/(2 - T]), 



(40) 



with 



2 - 77„ 



To obtain their strong coupling limit, we insert for 77^ and rj their perturbative expansions < pT\ ) and 
and reexpand the resulting ratios in power of ub up to the second order. This gives 



(41) 

81), respectively. 
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7 = l + 2(iV + 2)uB-4(7V + 2) 



The critical exponent v itself has the expansion 



- + {N + 2)uB-'i{N + 2) 



2(7V + 8) 



- (iV - 4) 



(42) 



2(iV + 8) 



(iV-3) 



(43) 



The strong-coupling hmits are, using p^/4 from Eq. (34) 



7 = 1 + 



(A^ + 2) 



2(iV + 8) - e(Ar - 4) 



_ 1 r (A^ + 2) 

^~ 2 \ ^ 2(Af + 8) - e(iV - 3) 



2 3(3iV+14) 
' (iV + 8)2 

2 3(3iV+ 14) 



(44) 
(45) 



(A^ + 8)2 

Their e-expansion are in agreement with _D = 4 — e results |^,^ • From these expressions we can recover 77 using the 
relation r/ = 2 — ^ /v. The result has now the correct e-expansion: 



iV + 2 2 

Tl=— 



(46) 



2(7V + 8)2 

This calculation of 77 via v and 7 was made in [^,^ to compensate the lack of a linear term in (|3^). Let us point 
out that, even if the e-expansion is not recovered, it is nevertheless hidden in a direct resummation of ( |38| ) to r/ = 0. 
To see this, we add a small dummy linear term ^u, to the defining equation (|3l|), leading to the expansion 

.4(A^-f-8)" 



V = Cub + 



8{N + 2)-C- 



(47) 



Using (|23D and (p4) , this leads to the strong-coupling value 



V 



P 



4 8(iV + 2) -4(7V + 8)C/e' 
which is zero for C = 0. Consider however the e-expansion of the right-hand-side performed at a finite C: 



(48) 



V 



4 4(iV + 8) 



1 



2(iV + 2) e 



(49) 



A^ + 8 C_ 

If we now take the limit C ^ 0, the right-hand-side starts directly like e^. Together with the lowest-order value 1 of 
p2/4, we obtain correctly (^6|). 



For consistency, the different two-loop results for 77, once from ( }44| ) and (|45| ), and once 77 = from (^8]) should not 
be too far from each other. This can indeed be verified by plotting the curves rj — 2 — j/v against a few values of N. 
The curves are all close to the 77 = axis for all N, approaching it for iV — > 00. 

Also for higher-loop orders, 77 could be obtained from the strong-coupling limit 7 and z/, or by taking a direct 
strong-coupling limit. Variational perturbation theory does not know which of these approaches should be better. 
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Ultimately, if we know enough term in the series expansion, the extrapolation to infinite order L should certainly 
become insensible to which function is resummed. 

One may wonder if it is possible to set up a unique optimal function of the critical exponents from which to derive 
the strong-coupling limit. The answer to this question would improve the theory considerably. 

Collecting the different results of this section, we have the D = 3 results 

1 



-1 + 2^/1 + 3(3iV + 14)/(iV + 8)2 
2N^ + 63iV2 + 540A^ + 1492 
(iV + 8)2(iV + 20) ' 
iV3 + 31iV2 + 2627V + 714 



(iV + 8)2(iV+ 19) 



2{N + 2) 
27V + 21 
2(iV + 2) 



1 



3(3A^+ 14) 



(A^ + 8)2 

(/V + 8)2 + 3(3iV+14) 



TV + 20 2(iV + 8)3 + 5{N + 8)2 + 3{N + 2)(3N + 14) ' 
1 3 3iV+14 



(50) 
(51) 
(52) 
(53) 
(54) 
(55) 



4:{N + 8) 4(iV + 8)3' 

where we also included the value of the renormalized coupling constant at the IR- fixed point. It is obtained from the 



one-loop series in u of the expansion (14) 



u — ub 



4(7V + 8) 



(56) 



We can restrict ourselves to one loop since it corresponds to a power u^. The two- loop calculation was however 



needed to get uo correctly, which itself enters 



With the help of (23), we obtain 



3 3iV + 14 



4 4(iV-f 8) 4(iV-f 8) 4(iV 



(57) 



Since only two critical exponents are independent |12 1^, all other can be derived from Eqs. (|50|)-(p4|). These two 
loop expressions are only a lowest approximation to the exact results. In the next section, we evaluate analytically 
the strong-coupling limit of the exponents at the three-loop level. 



D. Critical exponents from three-loop expansions 

The three-loop calculations are algebraically more involved. Moreover, as far as the critical exponents are concerned 
(we will see later that this is not necessarily true for the amplitude functions) the optimum of the function (|l^ ) is not 
given by the vanishing of the first derivative, but by a turning point, i.e., by the vanishing of the second derivative. At 
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the three-loop order, this implies that the parameter r in ( p^ is zero, leading to the three-loop strong-coupling limit 
result (p5|). It is this feature which renders the calculation analytically manageable, involving only a cubic equation 
for the determination of p (without r = 0, we would have had to solve an eight-order equation). In order to obtain ui 
to three loop, we generalize ( ^2[ ) to the same order, and find 

'{N + S,f 3(3iV + 14)" 



1 = -r, — log — - — ^ -UB + 16 

8(iV + 8)3 60(iV + 8)(37V + 14) 96C(3)(5iV + 22) + SSiV^ + 9227V + 2960 



(58) 

From this we extract the coefficients fi (i = 0, . . . , 3) of (kJ). The argument of the square root r turns then out to be 



negative, and the equation to be solved is (25). This is true not only for e = 1, but also for all e G [0, 1]. Since ([2£ 
is a cubic equation for p, there are three solutions, one of which is always negative, which we discard as unphysical, 
leaving us with two solutions. Only one of them is connected smoothly to the two-loop result. The purely algebraic 
form of the solution, generalization of the square root coming from solving (|3^), is somewhat too lengthily to be 
written down here. As a check, we have derived its epsilon expansion which reads 



2 



-e 



(59) 



3(37V+14) 96C(3)(57V-|-22)(iV + 8)-|-33iV3 + 2147V2 + 1264iV + 2512 
^ (TV + 8)2 ^ 4(7V + 8)4 

and leads to the correct e-expansion for lo — e/{p ~ 1): 

3(3A^+ 14) 2 96C(3)(5N + 22)(N + 8) + 3m^ + 5S8N^ + A288N + 9568 . , , 

"- = ^- (Ar + 8)2 ' + WTW ' ^ ^ 

which is the extension of ( |36| ) to the order e"^. This is to be compared with Eq. (17.15) of the textbook JT^ . 

The trigonometric representation is however compact enough to be written down here explicitly, at least for e = 1. 

Introducing an angle 9 and two coefficient oq, bo defined by 

[13776 + 47387V + N^{8N + 405) + 96(57V + 22)C(3)]^ 

= arccos ' 



( 



2 [106 + N{N + 28)] {{N + 8) [13776 + 47387V + 7V2(8iV + 405) + 96(57V + 22)C(3)]}^/^ 

1 

[2209664 + 1040160iV + 162982iV2 + 9683N^ + 184:N* + 672(7V -f 8)(57V + 22)C(3)]^/^ 
X |67181166592 + 640010403847V + 258933120007V2 + 56418284807V^ -I- 713027988iV^ -I- 54733044Af^ 

+ 27601577V^ + 883327V^ + 14407V^ 

- 192(7V + 8)(57V + 22) [4084864 + 19524807V + 3237067V2 + 200217V3 + 5147V^] C(3) 

+ 746496 [(7V -t- 8)(57V + 22)C(3)]^ |^ , (61) 



"° 446336 -I- 2132807V + 353347V2 + 21797V3 + 567V4- 864(7V + 8) (57V + 22)C(3) ^^"^^ 
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bo = 3y/{N + 8) [13776 + 47387V + N^{8N + 405) + 96{5N + 22)C(3)] 



X ^[2209664 + 10401607V + 1629827V2 + 96837V3 + 184iV4 + 672{N + 8){5N + 22)C(3)], 
the relevant root of can be written as 



P = -l + [106 + iV(iV + 25)]' - aobo cos (^^^ 



For the physicahy interesting cases = 0, . . . , 4, we obtain the values for D = 3 dimensions 



(63) 



iV 





1 


2 


3 


4 


p 


2.41829 


2.40384 


2.38683 


2.36910 


2.35157 




0.705073 


0.712332 


0.721069 


0.730405 


0.73988 


uj (Ref. |l|) 


0.812 


0.799 


0.789 


0.782 


0.774 



(64) 



Figure 1 illustrates the two- and three-loop critical exponent of the approach to scaling lo = el[p— 1) as a function 
of N calculated from (jsj) and (|6^, respectively. For comparison, we also give the three- loop unresummed result (|6^), 
evaluated at e = 1 and the theoretical values given in Tables 1 and 3 of [|l^. The latter are based on a five- loop 
analysis supplemented by a large loop order analysis. 



1 . 5 



. 5 



TV 



2 4 6 8 10 

FIG. 1. Two-loop (short-dashed) and three-loop (solid) critical exponent uj for different 0(A^)-symmetries. For comparison, 

the e-expansion (mixed-dashed) and the theoretical values of [Q (dots) are also given. 

Once oj is known to three loops, the other exponents and the strong-coupling limit u* of the renormalized coupling 
constant can be determined to the same order. To obtain u* , the two- loop expansion of u in powers ub is enough 
since it is of order 0{u^g). Recall that the three-loop expansion of u{ub) is needed only to calculate oj. From ( |l4|) we 



identify /i,/2,/3 and use (25) (since the argument of the corresponding r in ( |24| ) is negative) to obtain the critical 
value 
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e{N + 8)p{p+l){2p~l) 



2e{N- 



'(2p~l) 



12 [2(iV + 8)2 + 3e(3iV + 14)] 27 [2(iV + 8)2 + 3e{3N + 14)]' 
with p from (|64[). If we use instead the e-expansion of p given in (p9|), we obtain 

^ „ _ e , 3 3iV + 14^2^ 4544+ 1760iV+1107V2_337v3-96(A^ + 8)(5iV + 22)C(3)^3 



(65) 



4(iV + 8) 4(A^ + 8)3" ' 32(7V 

In the same way, we find the strong-couphng hmit of the critical exponents 7 and i>, as defined in ( ^ ) and 
together with ( pO| ) and (|3|), the latter two exponents being obtained from the the mass (|l|) and wave function 
renormalization, respectively. The three-loop perturbative expansions are 

"2(iV + 8) 



(66) 
(0) 



-f = l + 2iN + 2)uB-^{N + 2) 



- (iV - 4) 



+ 4(iV + 2) 



8(iV + 8)2 4(2iV2 - iV - 106) 



+ 194 + 7V(27V + 17) 



iy=- + {N + 2)uB-2{N + 2) 



2{N + 8) 



(TV -3) 



+ 4(iV + 2) 



4(iV + 8)2 2(27V2 + iV- 90) 



+ 95 + iV(A^ + 9) 



(67) 



(68) 



from which it is immediate to identify the expansion coefficients /o, . . . , /a which enter (p5[), to obtain 

e{N + 2) [e{N - 4) - 2(jV + 8)] p{p + l){2p - 1) 
^ 3 [8(iV + 8)2 - 4e(2Ar2 N - 106) + e2(2iV2 + 17N + 194)] 

8e(A^ + 2) [e{N - 4) - 2(iV + 8)]^ (2p - 1)^ 

27 [S{N + 8)2 - 4e(27V2 _ TV _ 106) + e2(2iV2 _^ ijjy ^ i94)]2 ' 
_ 1 _ e{N + 2) [£(jV - 3) - 2{N + 8)] p(p + l)(2p - 1) 
" ^ 2 12 [4:{N + 8)2 - 2e(2iV2 + iV - 90) + e2(iV2 + QN + 95)] 

e{N + 2) [e{N - 3) - 2(A/- + 8)]^ (2p - 1)^ 



(69) 



27 [A{N + 8)2 - 2e(2iV2 + ^ _ 90) + g2(^2 + g^y + 95)] 
where p for e = 1 can be obtained from (p4). The associated e-expansions can be obtained using (p9f). They read 



(70) 



7=1 



iV + 2 (A^ + 2)(7V2 + 22iV + 52) 



2(iV + 8) 



4(7V + 8)3 



{N + 2) [3104 + 24967V + 664iV2 _^ ^^^3 + - 48(A^ + 8)(5A'' + 22)C(3)] 3 



1 



iV + 2 

V ^ - -\ : r 

2 4(iV + 8) 



8(7V + 8)5 
(jV + 2)(Ar + 3)(iV + 2G) 2 
8(iV + 8)3 ^ 



{N + 2) [8640 + 59047V + 14127V2 + 897V3 + 2Ar4 - 96(Ar + %){hN + 22)C(3)] 



(71) 



(72) 



32(Af + 8)5 

Figures 2 and 3 illustrate the two- and three-loop critical exponents 7 and v, respectively, as a function of N . 
They are given by Eqs. (|l]) and (|69| ) for 7 and by Eqs. (^2|) and (|70| ) for j/. For completeness, we also plot the e- 
expansion ( [ti]) and (^2|) of the exponents, as well as the theoretical values quoted in Tables 1 and 3 of ||l^. Contrary 
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to the case of the critical exponent w, we see that the two- and three-loop critical exponents are very close together. 
This is a virtue of working self-consistently with uj obtained at the same loop order. In 1^,^, the extrapolated uj to 
infinite loop order was used instead. This implies that each loop order result for 7 and v was not very close to its 
asymptotic limit (contrary to what we get here). However, the extrapolation formula ( p^ ) works precisely for this 
case, and very precise extrapolated results for 7 and v could be obtained. In our present work, the critical exponents 
are not very far from their asymptotic limit, already at the two- and three-loop level. However, the extrapolation 
formula ( p^ ) cannot be used. It is not yet clear to the authors how it will be possible to extrapolate the five-loop 
results obtained using the present formalism. This question is left aside for a future work. We also note in passing 
that the e-expansion result is not too far from the values obtained in the strong-coupling limit. 



2 

1 . 75 
1 . 5 
1.25^ 

1 

. 75 
. 5 
.25 



N 



2 4 6 8 10 

FIG. 2. Two-loop (short-dashed) and three-loop (solid) critical exponent 7. For comparison, the e-expansion (short- and 
long-dashed) and the theoretical values of (dots) are also given. 



Ir 

0.8 
. 6<i 

0.4- 
0.2 



TV 



2 4 6 8 10 

FIG. 3. Two-loop (short-dashed) and three-loop (solid) critical exponent v. For comparison, the e-expansion (short- and 
long-dashed) and the theoretical values of (dots) are also given. 
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The critical exponent rj is obtained using 2 — j/v, with 7 and ly from ( |69[ ) and ([70|), respectively. It has the 
e-expansion 

(N^9\ ^ (N ^9\(N'^ - ^f\N ^979^ „ 

(73) 



^_ (TV + 2) ^2 (iV + 2)(iV2_ _ 272)^3 



2(7V + 8)2 8(iV + 8)4 

Let us also calculate directly the strong coupling limit of r] from its definition (| 



T] = 8{N + 2)u% - 8{N + 2){N + . 



+ l]u' 



(74) 



At the two-loop level, the result was zero. At the three-loop level, the calculation is different from that of 7 and f 
because there is no linear term in ub- This has already been discussed after Eq. (p6|): although we are working at 
the three-loop level, the optimum of the variational perturbation theory is not governed by a turning point but by an 
extrcmum for which the sign of the root r is the opposite to the usual case. The solution corresponds to r = — 1 and 
the optimum is ub = ^2/2/(8/3), so that 



77 



4 /I _ 32(2p-l)3(iV + 2) 2 



27 /I 27 [N 



(75) 



With (p9[), this leads again to the correct e-expansion (73). The difference between r] = 2 — j /ly and ( [75|) at e = 1 is 
illustrated on Figure 4 which also shows the direct evaluation of the e-expansion series (|7^) as well as the theoretical 
values quoted in Tables 2 and 3 of pSt]. 



.05 
. 04 
. 03+ 
. 02 
.01 

— . — . — . — , — . — . — . — , — . — . — . — , — . — . — . — , — . — . — . — ^ 
2 4 6 8 10 

FIG. 4. Two- loop (short-dashed) and three-loop (solid) critical exponent rj from the definition 2 — j/v. For comparison, 
the e-expansion (short- and long-dashed), 77 from the strong-coupfing limit of the direct (medium-dashed) series ( [tsI ) and the 
theoretical values of (dots) are also given. 

It is amusing to see that the e-expansion is the best approximation, followed by the strong-coupling limit of the 
direct series ([TSf). Comparing the different results, we see that they differ by about 30%. This is due to the absolute 
smallness of The error is small compared to unity. 
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To end up this section we also give the critical exponent 77™. Up to three loops, the bare perturbation expansion 
reads, from ( p^ ) and (|o|), 

"4(iV + 8)2 2(19iV+122) 



77,„ = 4(iV + 2)zlB-8(iV + 2) 



2(iV + 8) ^ 

— +5 

e 



u% + l&{N + 2) 



-3(5iV + 37) 



u%, (76) 



from which we deduce the strong-coupling limit with p from ( p^ 

e(A^ + 2)(2iV+ 16 + 5e)p(/9+l) (2/9-1) 4e(A^ + 2)(27V + 16 + 5e)3(2p - 1)^ 



3 [4(iV + 8)2 + e(387V + 244) + 3e2(5iV + 37)] 27 [A{N + 8)2 + e(38iV + 244) + ie^ibN + 37)] 



2 ■ 



(77) 



Its e-expansion is 



N + 2 {N + 2){im + 44) 2 {N + 2) [5312 + 26727V + 452iV2 - - 96(iV + 8){5N + 22)C(3)] 
N + S^^ 2(iV + 8)3 ^ ^ 8(7V + 8)5 



(78) 



The result (|7^) is analytically different but numerically close to that obtained via the scaling relation (pl|), implying 
rjm = 2 — j/^^, as illustrated in Figure 5. For completeness, the figure also shows the e-expansion (fz^ ) and the 
theoretical values quoted in Tables 2 and 3 of flq] . 




- — ^ — ^ — ' — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ' — ^ — ^ — ^ — ^ 
2 4 6 8 10 

FIG. 5. Two-loop (short-dashed) and three-loop (solid) critical exponent rim from the definition 2 — u~^. For comparison, 

the e-expansion (short- and long-dashed), rjm from the strong-coupling limit of the direct (medium-dashed) two- loop (js^) and 

three-loop (long-dashed) series ( [t?! ) and the theoretical values of Q (dots) are also given. 

We see a better agreement with the theoretical values quoted from when the exponent is evaluated in the 
strong-coupling limit of the direct series (39) and ([77|). This was also the same for the exponent 77. 

Collecting the different results of this section, we have the analytical form of the D = 3-dimensions critical exponents 
in the three-loop order 
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7 = 1 



{N + 2){N + 20)p(p + l)(2p - 1) 8{N + 2)(iV + 20)3(2p - l)^ 



3(2A^2 1497V + 1130) 



27(27V2 + 149iV+ 1130)2 ' 



_ 1 (iV + 2)(iV+19)p(p+l)(2p-l) _ (jV + 2)(jV + 19)3(2p- 1) 
^~2^ 12(Ar2 ^ 71^ _^ 531) 27(A^2 _^ 717V _,_ 531)2 

_{N + 2){2N + 21)p{p + l){2p - 1) 4{N + 2){2N + 21)3(2p - 1)^ 



3(47V2 _^ ii7iv + 611) 
32 (2p- l)3(iV + 2) 
2187 {N + 8)2 ' 
(iV + 8)p(p+l)(2p-l) 



27(47V2 + 1177V + 611)2 
2(7V + 8)3(2p- 1)3 



(80) 
(81) 
(82) 
(83) 
(84) 



12 [2{N + 8)2 + 3{3N + 14)] 27 [2{N + 8)2 + 3{3N + U)f ' 
where p is given in Eq. (|6^). For 77 and Tym, we took the strong-coupHng hmit of the direct expansions: Eq. ( [75[ ) for rj 
and Eq. (fz^ ) for Ty-m. These results have to be compared with the two- loop ones given in (|50|)-(|5^). 
For completeness, we give below the table of the critical exponents to three loops: 



TV 





1 


2 


3 


4 


7 


1.16455 


1.2338 


1.29426 


1.34697 


1.39307 


7 (Ref. II) 


1.1596 


1.2396 


1.3169 


1.3895 


1.456 


V 


0.587376 


0.623381 


0.654552 


0.681561 


0.705071 


V (Ref. @) 


0.5882 


0.6304 


0.6703 


0.7073 


0.741 


?7m 


0.311607 


0.421796 


0.509799 


0.580684 


0.638337 




0.0258218 


0.029917 


0.031452 


0.0315846 


0.03096 



They cannot compete with the five- loop calculation of l^~[^, [r^ , [l^ . However, our results are analytical, and already 
close to the asymptotic limit although we made no assumption about the large order behaviour of the theory. We 
consider this as promising. In a subsequent publication, we will present a numerical calculation up to five loops, with 
large-order behavior information included, of our self-consistent formalism. 



IV. CALCULATION OF AMPLITUDE FUNCTIONS AND RATIOS 

From now on, we shall focus entirely upon the D = 3-dimcnsions model. As we mentioned in the introduction, it 
is only for the critical exponents that the minimal subtraction scheme leads to the same resummed values both for 
Z? = 3 and D — \ — t. For this reason, it made sense to study the e-expansions of the critical exponents, which was also 
useful for comparing with calculations in 4 — e dimensions. The reason for this equality is the mass independence of 
the renormalization constants in this MS scheme. The mass independence implies a decomposition of the correlation 
functions into amplitude functions and power parts, for which the latter can be evaluated in the symmetric phase. 
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The amplitude functions, however, depend on being in the ordered or disordered phase. Moreover, the situation is 
complicated for > 1 by the presence of Goldstone singularities, most of which have to cancel at the end of the 
calculations: only the physical singularities, for example those occurring in the transverse susceptibilities, should stay 
at the end of the calculations. 

For this reason, apart from the three- loop work ]39| ], no three- or higher- loop calculation has been done for iV > 1 
below Tc, even numerically. The only relatively easy case is iV = 1 for which extensive numerical work has been done 
below Tc up to five-loop order plp^p^ . Above Tc, all N can be treated in the same way p0| , p^ , p9[ . In the latter 
reference the critical exponents 77 and rjyn have even been obtained to seven loops, with resummation performed in 

We have explained in detail in the first part of this paper that it is unnecessary to go to the renormalized theory since 
all results can be obtained from the strong-coupling limit of the bare theory. In the literature, the effective potential 
is given in terms of the renormalized quantities 1 40 39|j4^ . To apply our theory, we shall rewrite the expressions back 



in the bare form, using (fi^ 



A. Available expansions 

Let us list the most important available amplitude functions derived from the minimally renormalized model at 
Z? = 3 at vanishing external magnetic field hs- Up to two loops, they can be found in Rcf. Q]: 

• the square of the order parameter = {'i>%) below T^- 



f<t> — 



327rw 



— (160 - 82iV) + -{N - 1) ln3 
2"tt tt 



(85) 



the stiffness of phase fluctuations below Tc (some authors call this helicity modulus [^) T: 



/t — 



1 1 

8^ ^ 3 



54 



(2378 - 683iV) + 8(iV - 3) In 3 



(86) 



• the q'^ part of the transverse susceptibility xt- 



47V- 128 In 3 



(87) 



• the specific heat C above and below Tc'. 
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F+ = -TV - 2iV(iV + 2)u, 

F_^^-4 + 8(10 - N)u, 
2u 



(89) 



the isotropic susceptibility above Tc |52 



go 

f^^^l--{N + 2)u 



(90) 



• the amphtude function of the susceptibihty below Tc, which we obtain taking the inverse of the two loop 
numerical expansion given (up to five loops) in [^sf : 



/x_ = 1 + 18u + 164.44m2. 



(91) 



The latter quantity is restricted to = 1. 

From the series expansion of f^^ and /x , one sees that the relation 



(92) 



is satisfied to two loops. This is not a surprise: the bare helicity modulus, defined as T = 2dTB/dq^\q=o where is 
the free energy, can be shown (at least to two loops |4^) to be identical to M]^{dxT^ / dq^)\q=Q- This is a consequence 
of a Ward identity for the broken 0(A^)-symmetry below Tc- 

In Ref. po| , the perturbation expansions of the amplitude functions for the order parameter and for the specific 
heat have been carried to three loops. The additional terms are (we use the notation fj — ^ - /j'^w*): 



f(3) 



1 



2500A^2 + 65104Ar + 29056 + 8640(57V + 22)C(3) + 58320ci - lhTT^{im^ + 643iV + 499) 



10807r ^ 

- 180(647V2 6407V + 457) Li 



2 1--)- 80(194Ar2 16167V - 1675) ln3 + 16(8607V2 + 8357iV - 7867) ln2 



270(7V- 1) 



'8c2 + 32Li2 ( -i j +42 Lis ) - 64Li2(-2) + 21(ln3)2 + 16(ln2)2 - 96(ln2)(ln3) 



) = -m{N + 2) ( - — + 41n ^ 
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(93) 
(94) 



1 2 

= (10807V2 + 3464iV + 31120) - 128(57V + 22)C(3) - 864ci + -T:^{m^ + iV + 17) 

2T 3 

+ 216Li2 f -- j -32(47V+17)ln3 + y(31iV + 95)ln2 



+ 4(7V- 1) 



-8c2 + 16Li2 ( -ij +6Li2 Q ) -32Li2(-2) + 3(ln3)2 + 8(ln2)2-48(ln2)(ln3) 



(95) 



where Li2(a:;) = X]^o /^"^ is the dilogarithmic function [S^ , and ci and C2 are two numerical constants given by a 
single variable integration over elementary functions [p3lE0| : 



Cl 



C2 



dx 



Ve - 2a;2 



, 3 ,3 + 2: a; 
In - + In h T 

4 2 + X 2 + x 



In 



3 + a; 



a; 2 + 2: 
In — ■ — 



2-2: 



0.021737576333, 



4V2 



V2 



dx 



VlT: 



In- 



ln(l + 2:) 



1 + a; 



0.973771427. 
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(96) 
(97) 



For completeness, we give in Appendix^ some hints on how to obtain these amphtude functions. For the details, see 
Rcfs. [|9|,||. Our own contribution concerns the susceptibilities above and below T^. Using the three-loop integrals 
available in the literature |^,^ , we have been able to calculate analytically the thee- loop extension of the amplitude 
of the isotropic susceptibility f^_^ : 



-21 + 127r^ + 1281n- + 144Li2 ( 



as well as the three- loop amplitude function of the susceptibility below Tc, for = 1: 



(98) 



/x_ = 1 + 18u 



1480 



1072- 1 1664ci + 37r2 + 10480 In - + 36 Lia - - 



(99) 



Our analytical two- loop coefficient 1480/9 agrees with the numerical coefficient given in (^). We shall comment on 
the three-loop one later. The details of the calculation are given in Appendix ^ and 

B. Amplitude ratios 



Besides the amplitude functions, we shall also evaluate three important ratios: the amplitude ratio of the heat 
capacity, the universal combination Rc, and the amplitude ratio of the susceptibilities for A^ = 1. For a review 
of amplitude ratios, see [Q. The relevant equations for their determination is given in Appendix One of the 
best measured amplitude ratios was mentioned in the introduction: it is the amplitude ratio of the specific heat of 
superfluid helium above and below Tc, corresponding to A = 2. It can however be defined for all A^ and, using our 
notation, can be written as | |3^ , ^ 

A+ _ fb+Y fAiyB* + aFl' 



A- 



(100) 



where a and v are critical exponents and B* is the vacuum renormalization group function associated with the 
additive renormalization constant of the vacuum, evaluated at the critical point. It is known to five loops in the 
minimal subtraction scheme pO|p[ ] and reads, up to three loops 



uB = ^u + m{N + 2)u^ 



(101) 
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The ratio is equal to |3^: = 2i^P^/[(3/2) — 2h'P^], where P+ is a polynomial in u, related to the scale 

above Tc- Its analytical derivation is given in Appendix ^ and reads, up to three-loops: 



P+ = l- 2{N + 2)u + AiN + 2)u2 



+ 27(^ + 2) 



3 /I 

-3(63iV + 572) + 24(iV + 8)Tr^ + 4(437V + 182) In - + 288(iV + 8) Liz f -- 



(102) 



The experimental test for the validity of the strong-coupling expansion is to match (IOC) with (|l|) for N — 2. We 
shall see in the next subsection if this can be done. 

The ratio Rc is defined by the universal combination of amplitudes Rc = A'^ / Alj where r"*" and Am are 
the leading amplitudes of the isotropic susceptibility above Tc and of the order parameter below Tc, respectively. This 
ratio has been written in Ref. as 

(2iyP*)2-2/3 4:iyB*+aFl 1 



Rc 



(103) 



(3/2 - 2i.P;)-2/3 16^ /;/*^ ■ 

All the quantities have been defined previously, but for f3 which may be taken from the hyperscaling relation (3 — 
v{D — 2 + rf) / 2 = v{l + rf) / 2 in D = i dimensions. However, our own calculation for Rc gives a correction to (|103|): 



Rc 



(2;/P_;)2-''('D-2) AvB*+aFl 1 



(2;.P;) 



* \2-2/3 



(3/2-2i/P;)-'^(-D-2) 



167r 



fp*^ {i/2-2^PX)-^P \b 



4i^P* -I- aFl 1 



167r 



(104) 



<t>-ix+ 



Since this disagrees with (103), we give our derivation of this result in Appendix ^. We have verified that the numerical 



values coming from (103) and (104) do agree within 1%. This is traced back to the small value of the exponent ?/. In 



the following we shall however consider (104). We hope than the analytical discrepancy between (103) and (104) will 
soon be resolved. 

The third ratio to be investigated is the amplitude ratio of the susceptibilities for = 1. Such a ratio can also be 
defined for the longitudinal susceptibilities for A^ > 1. This is a nontrivial task requiring an appropriate description 
psf due to Goldstone singularities and this will not be investigated here. Using the notation of |3^,^^, the amplitude 
ratio can be written as 



'I Y 



b+ 



1v 



(105) 



where the ratio has been defined below Eq. (100), and where the quantities are restricted to A^ = 1. 

The question arises now to calculate the amplitude functions and ratios. As for the case of the critical exponents, 
we shall proceed also by order, starting with two loops. 
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C. Amplitude functions from two- loop expansions 



In order to apply strong-coupling theory to the amplitude functions (|85|)-(|88|), we must reexpand them in powers of 
the bare coupling ub using ( p^ up to two loops. The strong-coupling limit is then given by the general expression (|2^), 
with p2/4 given by Eq. (Isl) at e = 1. 



We start considering f^. To deal with a Taylor series, as assumed in the general theory in Section III A, we consider 



— (160- 82iV) + -{N- l)ln3 
277r TT 



(106) 



This series is special because the linear term in ub, is absent: the optimal value (|2^) is therefore given by u*^ = 0, 
and the two-loop value of u/^ in the strong-coupling limit is the same as the lowest-order value, which is independent 
of iV: 



327r' 



(107) 



It is worth pointing out here the effect of the special choice for Ajj in (J^ . We mentioned there that this coefficient 
did not have any influence upon the critical exponent. This is because the factor Ajj can be absorbed in to give 
Ub, implying the same strong-coupling limit. However, amplitude functions are ^^-dependent. In particular, for 
uf^, the chosen value A3 — l/(47r) has made the linear term disappear. One sees that this choice correspond to an 
optimalization: the zero order, the one-loop and the two-loop optimum values coincide. One expect then that the 
third- loop order contributes only to a small deviation from it. This is indeed the case, as will be shown in the next 
section, and confirm previous expectations ]3^j39[ |. 

The same situation holds for the amplitude function (^) of the susceptibility above Tc- The linear term in u being 
absent, the optimal value to two loops is independent of N and is equal to 



(108) 



The strong-coupling limit of the amplitude function of the stiffness of phase fluctuations and of the g^-dependent 
part of the transverse susceptibility can also be easily determined. The bare expansion is obtained combining 
(|^) and dll) to two loops: 



„ 1 MB 
"■^^=8 + T 



^(1802 - 7557V) + 8{N - 3) In 3 



/x^ 1 + -UB - -(ll/V- 58 + 96 ln3)u|. 



(109) 
(110) 
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The corresponding optima are given by (E3|) with from which we obtain 



u*f* 



1 



6(iV2 + 257V + 106) 



^ 8 (iV + 8)2 [7557V - 1802 - {432N - 1296) In 3] ' 
16(A^2 _^ 25iV + 106) 



f* = 1 

■'XT 



3(Ar + 8)2[llAr- 58 + 961n3]' 



(111) 



(112) 



The rcsuh ( pTl| ) has a pole for N = 2(648 ln3 - 901)/(4321n3 - 755) « 1.349, indicating that the strong-coupUng 
result is unrchable. We expect the pole to be an artifact of the limitation to two loops which disappears at the three- 
loop level. Since /x is not known to three-loops, we can only give plausible arguments for this expectation, suggested 
by the calculation of u*{Fl — F^) up to three loops in Eq. (|l36|), where a similar pole arises at the two-loop level 



but disappears for three loops due to the interplay of the coefhcients of the loop expansion. The trouble with (111) 
derives from the fact that the term of order in (109) change sign for the mentioned value of iV w 1.349, and at 



the two-loop level nothing can compensate this. This is in contrast with critical exponents which were observed to be 



alternating series in powers of ub- The result (112) for /^^ is smooth for all positive N. A more reliable result for 



/y than the singular (111) can therefore be obtained by combining (107) with (pT3) via relation (p2[), leading to 



XT' 



(113) 



Note that for > 4, far away from the pole, the two results (111) and (|ll3| ) agree within 2%. 

It is worth pointing out that an evaluation of the renormalized expression (|87|) at the critical point u* given by (pE 



leads to a result compatible with (112) within less that 1%. This is due to the fact that higher-order correction to 
the zero order result /^^ — 1 are small for all N. This is in contrast to and /x, where two- loop corrections are 
important. 

We now turn to the amplitude functions F± which enter the heat capacity above and below Tc. At the two- loop 
level, they are given by Eqs. ( |88| ) and (p^), respectively. With the relation between the renormalized and bare coupling 
constant (p^) to two loops, we have the expansions 



uF+ = -Nub + ^N^N + 14)m^, 
uF^ = 1 - 4ub + 8(A^ + 26)u%. 



(114) 
(115) 



With the help of (ga) and (M), we obtain 



N N'^ + 25N+106 

11 F = 

+ 2 {N + 8)^{N +U)' 

1 ^N^ + 25N+m6 
u FZ = - - 2 



(A^ + 8)2(A^ + 26)' 



(116) 
(117) 
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In 1 39 40], uF+ was not a good candidate for Borel resummation because its u-expansion lacking alternating 



signs of its coefficients. This problem is absent in variational perturbation theory since the expansion (114) in term of 



the bare coupling constant ub does have alternating sign. The latter is then expected to lead to a reliable result (116). 
This will be confirmed by the three-loop result of the next section. 

To apply the usual Borel resummation at the level of the renormalized quantities, Refs. |^,^ wrote the amplitude 
ratio of the heat capacity as 



1 - a- 



(118) 



instead of (100), and resummed u{F^ — F^) and uF^, avoiding the direct resummation of m_F-|_. For comparison, we 
give below the optimal value of the difference u(F„ — F+). It is is determined from the expansions ( ^8|) and 
Using (|l^, it yields 



u{F_ ~F+)^- + {N- A)uB - 2{N^ + ION - 104)u|. 



(119) 



Its strong-coupling limit is, from ( |23| ) and (34): 



..p. _ 1 (jV-4)^(jV^ + 25iV+106) 

^ " +'2 2{N + 8)^{N^ + 10N -lOA)' 



(120) 



The latter expression diverge for a positive value of iV = —5 + \/129 w 6.358. Then, the difference m(F+ — F_) is 
not the good quantity for the strong-coupling limit at the two-loop level. We should rather evaluate uF+ and uF^ 



separately in the amplitude ratio (IOC), instead of using the equivalent expression (|llq ). We shall see in the next 
section that the pole of u*(Ff — _F^) is an artifact of the two- loop calculation. A similar conclusion was also obtained 



for the strong-coupling limit of /x, see Eqs. (Ill) and (113). For ^ 4 and for N 3> —5 + vT29, the two-loop 



expansion (119) is alternating, and we expect that the strong-coupling result (12C) is reliable. As an indication for 



this, we compare ( |120[ ) with the difference of the optimized u*F^ values given in Eqs. ( |116[ ) and ( |ll7| ) 

, ^ 1 N'^ + 25N+106 f 2 N 

U*AF1 ^ -75 ; 

^ 2 (Af + 8)2 V^ + 26 2(A + 14) 

In Figure ^, we compare the two curves ( |120| ) and (121). 



(121) 
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FIG. 6. Comparison between the two-loop strong-coupling limit of u*{Fl — F^) and u* AF^. 



As far as the amplitude ratio of the heat capacity (JlOOj), or (|11S| ), is concerned, we still need to determine the 



strong-coupling limit of the renormalization group function B{u) of the vacuum ( 101 ) and of the polynomial P+ 



defined in (102). Because there is no contribution of the two- loop order to (101), its strong-coupling limit is 

N 

u B = u —. 



(122) 



Since the optimal two- loop result is identical to the one- loop result, it is clear that we may expect the large order 
limit L ^ (X) to differ only little from N/2. This has been confirmed in the five-loop resummation performed in 
and will be also seen in our three-loop calculation in the next section. 



The polynomial P_|_ given in (102) is evaluated in the strong-coupling limit using the same lines. The starting point 
is the expansion in powers of the bare coupling constant given in Eq. ( F4 ) of Appendix ^. Its two-loop part combined 
with Eqs. (|3|) and (||) leads to 

{N + 2){N^ + 25N +106) 



^+ (TV + 8)2 (27V + 17) 

The last amplitude we shall calculate using strong-coupling theory is . From (|2 



two-loop part of (D3), we find, with Eqs. ( p3| ) and (34), 

o2 



1 + 9- 



4 4532 



211/103. 



(123) 

at = 1 and from the 
(124) 



Combining with the unit value (|lOq ) of /*^, we have a ratio fx'_/.fx'+ identical to (124). However, this ratio might 
as well be determined as the strong-coupling limit of its perturbative expansion, instead of evaluating independently 
the strong-coupling limit of the numerator and the denominator. The relevant equation is given in Appendix ^ Using 



the two-loop expansion of (D6), we have, with Eqs. (p3|) and (|34|), 
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fx_ 

fx+ 



3 X 18^ 
T 1420 



751 
355' 



(125) 



D. Amplitude functions from three-loop expansions 



Some of the amplitude functions have been obtained to three loops. We now turn to their strong-coupling limit. 
This is done by applying Eqs. (^J), and (|6^ ) to the different amplitude expansions. 

We start with the amplitude function of the square of the order parameter. Combining the two-loop expansion (|85| ) 
with the three-loop term /j^^"* (p3|), and using also the relation between the bare and renormalized coupling con- 



stant (|lj), we have the three-loop expansion 



327r 



— (160- 82iV) + -(TV - l)ln3 
27tt tt 



— (160 - 82N) + -{N - 1) ln3 
277r TT 



u%. (126) 



From this, we read off the expansion coefficients /o, /i, /2, fs entering Eqs. (24), ( |25| ) and (p4[). Since the linear term 
/i vanishes, we have to follow the development below Eq. (^6|), adapting it to the present case. This development was 
done assuming a series with alternating sign since the expansions of the critical exponents had this property. Here, 
this is no longer true. Consider once more the derivation of the strong-coupling limit following from the optimal value 
of / = /o + f2U% + /sUb- ^b(2./2 + 3/3UI3) = 0. Two solutions are possible: u*g = and iig = —2/2/(3/3). The 
latter was relevant for the critical exponent rj. This does not mean that the other solution has to be rejected. In fact, 
looking at the nature of the extremum (minimum or maximum) , we see directly that the first solution corresponds to 



du% 



2/2 



(127) 



ub—u%>—0 



while the other leads to 



dill 



-2/2. 



(128) 



iiB=UB= -2/2/ {3 fa) 



If one solution is a maximum, the other one is a minimum. Looking for the sign of /2 in Eq. (126), we see that it is 



positive for N < TV^, with — (80 — 27 In 3)/ (41 — 27 In 3) w 4.43992, corresponding to a maximum, and negative for 
N greater, corresponding to a minimum. Variational perturbation theory at loop order L > 1 says nothing about the 
nature of the extremum. It might be a minimum or a maximum. In quantum mechanics, this has been explained in 
the book j^]. In quantum field theory, the exponent ?/ illustrates this: we had chosen the maximum (recall Eq. ([zs])). 
In this way, the e-expansion was obtained. Taking the solution u*g = 0, corresponding to the minimum, we would 
have obtained the three- loop result 77 = 0. The lack of reproducing the e-expansion gives a hint that the maximum 
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solution has to be chosen. In the case of we can also argue that the maximum solution has to be chosen, although 
there is here no e-expansion available, by definition of the model. However, at the point were /a = 0, we have to 



recover an optimization problem of a quadratic equation in ub, see Eq. (126). We know for this function that, because 
no linear term is present, the strong-coupling limit is u* = l/(327r). This implies that Ug = Oa,t this point, i.e., 
that the maximum solution has to be chosen. By continuity, this remains true in a neighborhood. The nature of the 
solution can only be changed when both solutions are equal, i.e., for N smaller than its value making /2 vanish. 
Below N^, we can imagine that we have an interchange of solutions, and that the minimum has to be chosen. In this 
case, we would have u* = l/(327r) for all N. If we decide to keep the maximum for all N, which we could prove 
to be true only for N > N^, this would imply that u*g = —2/2/(3/3) has to be chosen below and u*^ — above. 
Below we have f* — fo + 4/|/(27/|), as was the case for the critical exponent r/, while above N^, the solution is 
/* = /o- The strong-coupling limit of /^ to three loops is then 

' |/j^'-8(iV-f 8)[(160- 82iV)/(277r) + 2(iV-l)(ln3)/^]| ^ ' ' 

with p given by Eq. (|6^), and where 8 (a;) is the step function of Hcaviside, being equal to 1 for a; > and being 
vanishing for x < 0. As mentioned, we cannot be assured that for N < the maximum has still to be chosen. The 
possibility that the three-loop result is identical to the two- loop result remains. Would the above analysis not be 
performed, i.e., choosing the minimum solution everywhere, we would have obtained (|129| ) without the step function, 
meaning the presence of a pole at the vanishing of the coefficient of the cubic term in ( |126| ), i.e., for N « 4.92915. We 
have checked that the solution is sharply peaked near this value so that it would appear that the optimal value is valid 
almost everywhere. In fact, since the pole gives a very peaked contribution, a calculation at fixed integer value of N 
would have missed it completely, making one to believe that the resummation was correct. But this would not be 
true, the true solution being ( |129| ) everywhere. We give in Figure [t] the comparison between our two- and three-loop 



results. Our values for N < lie above the two-loop result l/(327r) obtained in (107) 
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FIG. 7. Comparison between the two-loop (short-dashed) and three-loop (solid) amplitude function of the order parameter. 
The resummed values |^ obtained using a Borel resummation are indicated by the dots for values of A'' available. 

This is also the case for the resummed values given in [M for TV = 2, 3 as can be seen in the following table: 



N 





1 


2 


3 


4 


N > N4, = (80- 271n3)/(41 - 271n3) 


u*f; (2 loops) 


l/(327r) 


l/(327r) 


l/(327r) 


l/(327r) 


l/(327r) 


l/(327r) = 0.00994718 


u*f; (3 loops) 


0.0105523 


0.0102518 


0.0100884 


0.00999735 


0.00995195 


l/(327r) 


u*f; (Rcf. 0) 






0.010099 


0.00997 







The agreement between our two- and three-loop order, and between our work and |10| |, is excellent. It is due to the 
fact that the term of order zero contains almost all information on this amplitude. 

The three-loop amplitude functions uF+ and mF_ are given by Eqs. (|8^), (p^, ( p9[ ) and (|9^). As in the case of the 
previous amplitudes, the present expansions may not be alternating. This may make the argument of the parameter 
r in Eq. (^J) positive, so that ( p4|) has to be used to obtain the strong-coupling limit rather than (^5|). However, (|25| ) 
remains correct for all N for uF^, while the alternating property is lost for uF_ for N > 40. Since the physical cases 
corresponds to = 0, 1, 2, 3, 4, we can ignore the alternative ( p^ ) and ( |25| ) is used throughout. Using the relation (|lj) 



between the bare and renormalized coupling constant, the three-loop bare extension of Eqs. (114) and ( |115[ ) are 



uF+ = -Nub + 2N{N + 14)u| + \f^^ - 2AN{7N + 46) 



= i - 4uB + 8(iV + 26)u| + - 480(3iV + 22) 



(130) 

3- (131) 

This allows to identify the appropriate /o, /i, /2, /s functions to enter Eq. (|25|). In the strong-coupling limit, we obtain 

2N^{N + 14)p(p -f l)(2p - 1) 2 [2N{N + 14)]^ (2p - 1)^ 



u*Fl = 



6 



- 2AN{7N + 46) 27 \f^^^ ~ 2AN{7N -f 46) 



2 ' 



(132) 



u*F* 



1 32(iV + 26)p(p+l)(2p-l) 2 [8(iV + 26)]^ (2p - 1)3 



F^'> - 480(37V + 22) 



27 



- 480(3iV + 22) 
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(133) 



with p from Eq. (c 



Figures @ and Hshow the comparison between the two-loop (116) and (117) results of the previous section and the 



corresponding three- loop (|132|) and (133) results, as well as a comparison with values given in [p9[, when available. 





N 



FIG. 9. Comparison between the strong-coupling limit of the two-loop (shot-dashed) and three-loop amplitude func- 
tion (solid) u* Fl. The resummed values obtained using a Borel resummation are indicated by the dots for values of 
available. 



To be more precise concerning the comparison with |39|, we give the appropriate values of u*Fl in the next table: 
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N 





1 


2 


3 


4 


u*Fl (2 loops) 


0.372596 


0.379287 


0.385714 


0.391707 


0.397222 


u*Fl (3 loops) 


0.374166 


0.378474 


0.384065 


0.389883 


0.395484 


u*Fl (Ref. 11) 




0.3687 


0.384 


0.387 





From this table, we see that that the strong-coupling limit results for u*F^ at the three-loop level differ only a 
little from their two-loop counterpart. This was also see on Figure ^ and on Figure ^ for u*F!^. For = 1, we can 
also infer from the table that the results coming from variational perturbation theory and from a Borel resummation 
pgf are not in excellent agreement, not even within the error-bars of the latter: u*F*(N = 1) = 0.3687 ± 0.0040. The 
agreement is however recovered for the values N = 2, 3. 

For u*F^, there is no available comparison between our work and others. The authors of ||3^ could not performed a 
reliable Borel resummation, presumably because of the lack of an alternating series. A comparison is however possible 
for the difference u*(FZ — F^). We have seen in the previous section that the two-loop evaluation of this difference in 
the strong-coupling limit did not work well in our case because the second-order term in the bare expansion changes 
sign for some value of N. Let us see how the situation changes at the three-loop level, which has the expansion (see 



Eqs. dlSOj) and (|131| )) 



1 



- F+) = ^ + {N ~ 4)uB - 2(iV^ + 107V - 104)m^ + - -|- 24(77^" - UN - 440) 



Ug. 



(134) 



The coefficient of the second-order term vanishes for TV = —5 + vT29. This is not anymore a problem since there 
is a three- loop order term preventing a l//2-behavior, see the comparison of (p3|) and (Eq). The coefffcient of the 



three-loop order can itself vanish. For (134), this happens for N = N ^ 10.5324. Since (^4|) and ( |25D imply a behavior 
like I//3, it is legitimate to wonder about poles. The answer is simple: if the coefficient of the three-loop term 
vanishes, then the problem is formally equivalent to evaluating the strong-coupling limit of a two-loop series. The 
coefficient of the linear and quadratic terms are however different from the two-loop result since the linear term has 
a factor p{p + l)/2 instead of p and the quadratic one a coefficient (2p — 1) instead of a factor 1. We conclude that 
when the three-loop term vanishes, the strong-coupling limit should be well-behaved, giving a smooth curve around 
N. This discussion shows that the function r in ( p4[ ) is not always zero here because r contains the coefficient /2 
of the two-loop term, and its zero govern the behavior of the solution (|2^). We note here the important following 
point: the positive square root +r was chosen in ( p^ ) in order to match with a vanishing f^. We explained, and this 
was used when evaluating rj, that the negative root might play a role as well. For 77 to three-loop order, we had a 
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negative r. Let us see what happens for /2 = 0. The expansion to be optimized is / = /o + fius + fsu^, such that 
we have to solve /i + ^/^{ug)^ = 0. For the critical exponents, the signs of /i and /a are the same because the series 
are alternating. For this reason, this equation has no real solution, and wc must solve the turning-point equation 
Gf^Ug — 0, which is iig — 0, leading to the optimized result /* = /q. For the amplitude functions, we have already 
seen that the alternating property is not necessarily true, so that the solution — —fi/iifs) is real. At the point 
where /2 vanishes, we can see that the optimal value is 



(135) 



the positive or negative sign being chosen to get a continuity of the solution around /2. 

For the difference function u{F^ — it is possible to follow exactly the strong-coupling limit as a function of 
N. Depending on N, there are four different solutions: below iVi « 2.48527 and above « 16.6066, the argument 
of the square root of r is negative, and one uses Eq. (^5|). For N e [7V2,iV3[, one uses Eq. ( |2^ ) with the positive root 
(r = |r|), where N3 — —5 + \/129 w 6.35782 is the value of N for which /2 changes sign. The zero of /s lies within 
the same region. Finally, the last region is within the range N e]iVi,7V2], for which we use Eq. (^4|), but with the 
negative root r — — \r\. More precisely. 



u*{Fl -F*) = - + 



1 2(iV-4)(iV' 



6 



-|-107V-104)p(p+l)(2p-l) / _ 
'^^^ +24(7 -UN -UQ)] V 



p(y _ p(^^) _^ 24(77V2 _ i4jv _ 440) 
2[2(A^2 ^ iQ]^ _ i04)]3(2p - 1)3 



27 



p{3} p^o, ^ 24(77V2 _ UN - 440) 



p(3) 



(136) 



with r = for < A^i and iV > A3 , and r the negative or positive square root of 



(A^ - 4) f[^^ - Fl"' + 24(7A2 - UN - 440) pip + 1) 



n(3) 



r-" 1 - 3- 



(137) 



2[2(A^2 _^ iQN - I04)]2(2p - 1)2 
for A^ eliVi, A"2] or A' e [N2, N3I respectively. 

Thus, the pole in A^ of the two-loop approximation to u*{Fl — F^) was only an artifact of the low order. At the 
three-loop level, the singularity is avoided by the interplay between the different possible solutions of ( |2^ ) arising 
from the different branches of r: r = 0,±|r|, with \r\ to be identify with the hmction r defined below Eq. ([24|). This 
possibility was not exploited in previous works because of the alternating signs for the critical exponents. (See 
however rj which required r = — 1 for the strong-coupling limit of ([t^).) 

In Figure ^0[ we show the strong-coupling limit of u*{Fl — F^). For comparison, we also give the direct difference 



u*AF^ between u*F* and u*F^, as obtained from Eqs. (132) and (132), as well as its two-loop counterpart ( 121 ). 
The range for A^ has been increased to 30 in order to investigate the regions delimited by iVi, N2 and A3. 
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N 

5 10 15 20 25 30 
FIG. 10. Comparison between the strong-coupling limit of the two-loop (shot-dashed) and three-loop amplitude func- 
tion (solid) u*AF± — u*Fl — u*F^. The three-loop (long-dashed) evaluation of u*{Fl — F^) is in better agreement with 
the resummed values (dots) obtained in pOl using a five-loop (A^ — 1) or three-loop (TV 7^ 0) Borel resummation. 



For the direct difference, the changes brought about by the three-loop is very small, as before in Figs. || and ^. The 
difference between u*{Fi^ — FJ) and u*AF^ is however somewhat larger. To facihtate the comparison, the following 
table should be of help: 



N 





1 


2 


3 


4 


M*AF| (2 loops) 


0.372596 


0.433608 


0.485714 


0.530258 


0.568519 


u*AFl (3 loops) 


0.374166 


0.432926 


0.484899 


0.530224 


0.569615 


u*{Fl - F*) (3 loops) 


0.374166 


0.421864 


0.461436 


0.489995 


1/2 


u*{Fl-Fl) (Ref. Hi) 




0.4179 


0.461 


0.498 





The simple value 1/2 for the case = 4 comes from Eqs. (136) and (137): for = 4, r is vanishing, meaning the 
third term of ( |136| ) does not contribute. Since the second term is also proportional to — 4, only the zero loop order 
survives for the Higgs case. Our results for u*{Fl — F^) are in good agreement with the Borel results of Ref. [ p9[ . 
This is probably not a coincidence since we now resum the same function as they did. We note however that, for the 
Ising model {N — 1), we are not within the error bars of We have already noted this for the strong-coupling 
limit of u*FZ. 

Before cloturing the investigation of u(F_ — -F+), we recall the case of fr, whose direct two-loop strong-coupling 



limit gave Eq. (Ill), exhibiting a pole. We know the strong-coupling limit should not have been far from ftj./S, see 



the discussion leading to Eq. (113). We have shown in this section how a pole in u*{Fl — F^) at the two-loop level 
might disappear at the three- loop one. This is probably the case for fr- It would be very useful to get its three- loop 
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order. 



We can now turn to the strong-coupling limit of the renormalization group constant of the vacuum B{u). Its 



three-loop value has been given in Eq. (101). 

Upon inserting the relation between the renormalized and the bare coupling constant dlj), we obtain 



N _ 



uB{u) = ^—UB - 2N{N + 8)ul + N{8N^ + 167N + 686)u|j. 



(138) 



The series is alternating and behaves as for the critical exponents. No subtleties arises here as in the case of and 
u*{Fl — F^). In particular, the argument of the square root of r in Eq. (|4|) is negative for all N and we have to 
work with ( ^ ) . Using ( p5|) , the strong-coupling limit is 

N{N + 8)p{p + l)(2p - 1) 16iV(iV -t- 8)3(2p - 1)^ 



u B — 



(139) 



6(8iV2 ^ i67iV -I- 686) 27{8N'^ -f 1677V + 686)2 ■ 

This result is plotted in Figure |ll| together with the two- loop result u*N/2, see Eq. (122). We also indicate the 
approximate result u*B* — u* N/2 with u* from the three-loop expansion (Isl). There is no visible difference between 



the latter and (139) 




N 

2 4 6 8 10 

FIG. 11. Comparison between the strong-coupling limit of the two-loop (short-dashed) and three-loop (solid) renormalization 

group constant of the vacuum u* B* . For completeness, we also give the approximate three-loop (long-dashed) result u*^-^N/2. 

A five-loop calculation using Borel resummation is included (dots) for values of A'' available. 

To facilitate the comparison between the different approximations, we recapitulate the numerical results in the next 
table: 
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N 





1 


2 


3 


4 


u*B* (2 loops) 





0.0226337 


0.04 


0.0535312 


0.0642361 


u*B* (3 loops) 





0.0221074 


0.0391089 


0.0523643 


0.0628447 


u*B* = U(3)A^/2 





0.0219975 


0.0388885 


0.0520441 


0.0624386 


u*B* (Ref. [||) 





0.020297 


0.0363919 


0.049312 





For the comparison with [p9|, we have multiplied their five- loop results for B* with their five-loop u* . These five- loop 
results are within 7% from our three- loop strong-coupling calculation. This confirms that B* w N/2 to all orders. 
We shall however see in the next section that this 7% difference leads to a non-negligible difference in the universal 
combination Re- 

The next quantity we shall resum to three loops is the polynomial P+. The three- loop bare expansion of the 



renormalized P+ was given in Eq. ( |F4| ) and resummed to two loops in (123). The series in the bare coupling constant 
is alternating, and behaves as for the critical exponents. The strong-coupling limit of to three loops is then given 

by 

9 (iV + 2)(2iV+17)p(p+l)(2p-l) 



P_ 



1 + --. 



+ ^ ' 2 [-3(367V2 + 837iV -I- 3920) + 24(7V + 8)7r2 + 4(437V + 182) ln(3/4) -f 288(7V + 8) Li2(-l/3)] 

[N + 2){2N + nf{2p-lf 



54- 



(140) 



-3(367V2 -I- 837iV + 3920) + 2A{N + 8)7r2 + 4(437V + 182) ln(3/4) + 288(iV -I- 8) Lia (-1/3)]'' 
with p from Eq. (|6^. 

In Figure |l^, we compare ( |140D with the two-loop result from Eq. ( |123| ) . Almost no difference is found between our 
two- and three-loop expansions. 
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FIG. 12. Comparison between the strong-coupling limit of the two-loop (short-dashed) and three-loop (solid) polynomial 
P+. The values pol obtained using a five- loop Borel resummation are indicated by the dots for values of A'' available. 
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For a better comparison, we quote the numerical values for TV = 0, 1, 2, 3, 4 in the next table: 



N 





1 


2 


3 


4 


PI (2 loops) 


0.805147 


0.74269 


0.695238 


0.658642 


0.63 


PX (3 loops) 


0.807683 


0.745874 


0.698901 


0.662717 


0.63447 


P; (Ref. |§) 




0.7568 


0.7091 


0.6709 





The two- and three-loop results agree within 1%. The results agree fairly well with the five- loop Borel resummation 
performed in j3^] . We shall however see later that amplitude ratios depends crucially on the exact value of P'^ . For 
this reason, our three-loop calculation is probably not precise enough. We shall present in the last section a numerical 
five-loop strong-coupling evaluation of to more firmly settle this statement. 

To conclude this section, we discuss the amplitude of the susceptibilities above and below Tc to three loops. We 
already know from the previous section that the two-loop amplitude above Tc is identical to the order zero: /*^ — 1, 



see Eq. (108). As for the case of uf^, we then expect a very small deviation from the zero-order value as well as a very 



smooth A^-dependence. The series to evaluate in the strong-coupling limit is given in Eq. (C9). It is alternating and 
behaves like the series of the critical exponents. Moreover, with a vanishing linear term, but with a negative coefficient 
of the quadratic term, the solution of the optimalization problem is at variance with the case of the exponent 77 or 



the amplitude if, as for these quantities, we admit that the solution is a maximum. From Eq. (127), we determine 
that the optimal value is u*^ ~ 0, so that 



f* = 1 



(141) 



remains true at the three-loop level: The amplitude of the susceptibility above Tc at the three-loop level does not 
depend on N\ This is in contrast to [5^ , where a A^-dependent fit, using Borel resummation, has been performed. 
Because our resummed value up to three loops is f^^ = 1, it is tempting to conjecture that this is true for all orders. 
However, contrary to the case of 77 and J^, we have here no argument to tell that the maximum has to be chosen 
instead of the minimum. Only when going to higher orders, then having more expansion coefficients, can we decide 
which solution is the right one. For this reason, wc also mention below this other solution, which differs from unity 
for at most 2.5%: 

48668(7V-f-2)3(2p- 1)3 



f* = 1 



[27(A^ + 2)(A^ + 8)]2 [-113-h 127r2 + 1281n(3/4) + 144Li2(-l/3)]^ ' 
The comparison between the two curves is given in Figure 03, as well as a comparison with the fit 



(142) 
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(143) 



taken from Table 1 of [|3|, with = 9.68 {N = 1), 11.3 {N = 2), 12.9 {N = 3), combined with the five-loop u* of 
Ref. M. 



1 . Olr 



0.99- 



J?* 



x+ 




. 98 



2 4 6 8 10 ^ 

FIG. 13. Comparison between the two-loop strong-coupling limit (short-dashed) of the amplitude /^^ of the susceptibility 



above and the second possible solution (142) at the three-loop level (solid). The values p3| obtained using a five- loop Borel 



resummation (dots) are given for values of N available. 
More precisely, we have the table: 



N 





1 


2 


3 


4 


/x+ (2 loops) 


1 


1 


1 


1 


1 


f*^ (3 loops) from (141) 


1 


1 


1 


1 


1 


f*^ (3 loops) from (142) 


0.979543 


0.976298 


0.975082 


0.974977 


0.975472 


f*^ (Refs. ||,||) 




0.976791 


0.9748331 


0.9740978 





The fact that our three- loop calculation (142) agrees very well with Refs. [p9^j53l might be an indication that (|l42) 



should be preferable to (141). However, from a variational perturbation theory point of view, nothing can be said. 
Only the determination of the next order might resolve the ambiguity. 

Finally, we determine the strong-coupling limit of the amplitude of the susceptibility below for N — I. We have 
checked that the parameter r in (e3) is zero, i.e., we have to work with the turning-point equation (Eq). Applying it 



to (D4), we have 



4352p(p-hl)(2p-l) 



3 [19904 - 11664C1 -I- Sn"^ + 10480 ln(4/3) -I- 36 Li2(-l/3)] 

164852924416(2p- 1)^ 
19683 [19904 - 11664ci + Stt^ + 10480 ln(4/3) + 36 Li2(-l/3)]^ 



(144) 
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with p from Eq. (|6^. Numerically, this is evaluated as /*_ « 2.09387, to be compared with the two-loop (124) result 
211/103 w 2.048544. They agree within 3%. 



For the ratio (105), the calculation of fx-l fx+ needed. Its strong-coupling limit can be determined using the 
individual strong-coupling limit of the numerator and the denominator. In that case, the ambiguity on /*_ at the 



three-loop level is relevant. According to the choice /*^ = 1 from (141) or /*^ sa 0.976298 from (142) with iV = 1 
we have 

fx 



f* 

f* 

■'x+ 



2.09387, 



2.1447, 



(145) 
(146) 



respectively. 

The strong-coupling limit of the ratio fx-/ fx+ '^^^ ^^^^ be computed from its perturbative expansion. It has been 
derived in Appendix^ see Eq. (D6). The strong-coupling limit reads 



fx_ 

fx+ 



= 1 



1420p(p+l)(2p-l) 



3 [19184 - 11664ci -1- 997r2 + 94561n(4/3) -I- 1188 Li2(-l/3)] 
5726576000(2/9- 1)^ 



(147) 



729 [19184 - 11664C1 -I- 997r2 -|- 9456 ln(4/3) + 1188 Li2(-l/3)] 
Its numerical value is 2.11227, to be compared with the two-loop ( ]125| ) result 751/355 « 2.11549. The three-loop level 
is in very good agreement with the two-loop result, within less than 0.2%. However, this by no means signify that the 



asymptotic limit has been reached, and the ratios (14£) or (146) might be closer to the true ratio than (147). This 
is due to the fact the even and odd orders are on different converging lines because odd (even) terms come from an 
extremum (turning-point) condition or vice-versa. To see the speed of convergence, it would be necessary to compare 
the fourth-loop order with the two-loop order, and the fifth-loop order with the third-loop order. 



E. Amplitude ratios from two- and three-loop expansions 

We have now everything in hand in order to compute the ratio of the heat capacity /A^ , the universal combina- 
tion Rc and the ratio of the susceptibilities r+/r_. This section is restricted to a full too- and three-loop calculation. 
In order to improve the ratios, we shall break our rule of being selfconsistent in the next section and use there the 
maximum information available. 



We start with the heat capacity A'^ /A . Since we have a preference for (|100|) over (118), we shall work with the 



separate strong-coupling limit evaluation of u*Fi^ and u*F^. We have checked that the effect on the ratio A+/A- is 
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negligible. The exponent a entering it is calculated from the two- or three- loop result for i/ given in (^2|) and (|8l|), 
respectively, using the hyperscaling relation a = 2 — Dh'. 
Combining the different results derived previously, wc have 



N 





1 


2 


3 


4 


A+/A^ (2 loops) 





0.489106 


0.843065 


1.12691 


1.37015 


A+/A_ (3 loops) 





0.491088 


0.862439 


1.16719 


1.43243 



Regarding the fact that the critical exponent a is far away from its asymptotic limit (it is still positive for N — 2, 
while the shuttle experiment shows clearly a negative value), the results of this table are promising: For N — 2, 
we obtain Aj^/A^ w 0.862439 at the three-loop level, while the shuttle experiment ||2^ gives Aj^/A^ w 1.0442, see 
Eq. (^. We shall see in the next section that working with asymptotic critical exponents leads to a better agreement 
with experiments. 

The next ratio we examine is (|l04| ), the universal combination Re- The results are best displayed in a table: 



N 





1 


2 


3 


4 


Rc (2 loops) 





0.062474 


0.124819 


0.184355 


0.239967 


Rc (3 loops, = 1) 





0.05944 


0.121628 


0.182413 


0.239691 


Rc (3 loops, from (142)) 





0.060883 


0.124736 


0.187094 


0.245718 



We see an overall agreement between the two- and three-loop results. We have also checked that the ratio Rc 
calculated with the formula ( |103| ) used in |^] is within less than 1%. Moreover, our results are in agreement with the 
values Rc{N — 2) = 0.123 and Rc{N = 3) = 0.189 given in Ref. Since we expect that using the true critical 

exponent leads to a better ratio A^/A^^ it is important to see how Rc evolves. Will the agreement with be lost? 
This issue is investigated in the next section. 

To end this section, we study the ratio of the susceptibilities r^/F- for the Ising model, see Eq. ( |105| ). 

The two-loop result for the amplitude ratio ( |l05|) is, with v — 56/90 from (|5^), with = 127/171 from (|l23| ) and 
with = 1: 



F_ 



211 



1v 



211 /14224\ 



56/45 



3.69171. 



(148) 



103 V 3 - ^vPl j 103 V 8861 j 
This is still far from the value w 4.7 quoted in the literature ||l^,p3|. A small improvement is obtained using the direct 



strong-couphng evaluation of fx-l fx+ of Eq. (|125| ): 
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This value is still far from the expected ratio 4.7. However, the ratio depends sensibly on the value of the critical 



exponent v. For example, using v = 0.63, we increase (149) to T^/V^ = 4.002. The sensibility is also seen when 
calculating the three-loop value of the ratio: 



3.88785, 



(150) 



where the ratio [fx- 1 fx+)* been obtained from (147). 



F. Amplitude ratios using maximum information 



Up to now, we have followed the strategy to make a fully consistent two- and three-loop calculation. The comparison 
between the two- and three-loop amplitude functions has made us believed that the resummcd values are close to 
the extrapolated limit L — s- oo, although one has to take care that odd and even approximations are on different 
converging lines. For the critical exponents, it is primordial going to the asymptotic limit. For example, we have 
a{N = 2) still positive at the three-loop level, while the shuttle experiment, see second reference of |^ and Eq. (|^), 
shows a value of a(iV 2) = -0.01056. 

In this section, we shall relax our constrain of working only with two- and three-loop quantities and will take the 
maximum available information, i.e., our three- loop result for the amplitudes and extrapolated, or experimental, value 
for the critical exponents. We shall also see the effect of using uB to five loops. 

Except for a(N = 2) that we took from the shuttle experiment the exponents are taken from the D — "i tables 
of lH^, i.e., we are working with, for = 0, 1, 2, 3, 4: 



V = 0.5882 0.6304 0.6703 0.7073 0.741 
a = 0.235 0.109 -0.01056 -0.122 -0.223 



(151) 



Combining the three- loop strong-coupling limit of the amplitudes performed in section IV D with these exponents, we 
obtain, for A^jA^: 



N 





1 


2 


3 


4 


A+/A^ 





0.543406 


1.04516 


1.54386 


2.0444 


A+/A^ (Ref. H) 





0.540 


1.056 


1.51 





We have checked that the increase from the three-loop value (for N = 2, this ratio was 0.862439) is mainly due 
to using the correct a. For example, with the correct a but still using the three-loop of (|8l|), we would have 
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obtained, for = 2, a ratio 1.04711. It also does not depend too sensitively on using the five-loop strong-coupling 
limit of u*B* and P^, neither on using u*{Fl — F^) instead of the separate calculation of u*Fl and u*F^. For 
example, playing with all these quantities, the ratio, for N — 2, could be changed from A+/A^ = 1.04516 to, at most, 
A^/A^ ~ 1.049, depending which quantities are taken to five loops. A complete numerical study of this ratio, using 
variational perturbation theory up to five loops, will be presented elsewhere p6[ . 

For N ^ 2, our result 1.04516 coincide remarkably well with the shuttle experiment, see second reference of [p8|. 



For iV = 1, we have 0.543406, which agrees reasonably well with the values quoted in Table 5 of |18|, values which 
are both experimental and theoretical. In Ref. js^], the authors obtained 1.056 for N = 2. Their Table 4 make a 
comparison between their result and other works and experiments, for A'^ = 1, 2, 3. We see that the agreement is good. 
In our table, we have listed only the values calculated in the work since the model is the same. 



For the universal combination Rc, we obtain, using the five- loop critical exponents (151) and the three-loop 



amplitudes of section [V D 



Rc = 0, 0.0616257, 0.130341, 0.201404, 0.270882 (152) 

for A^ = 0, 1,2,3,4. 

Here also, we have checked that the main effect is due to choosing the correct a. Working with v at the three-loop 
level only modifies slightly the result. While working with the true exponents for the ratio A+/A^ had considerably 
improved it, making it coincide with the experimental values, we see for Rc that the values of the previous section, 
with a wrong a were in better agreement with the quoted values in |^^: Rc — 0.123, 0.189 for iV = 2, 3, respectively. 
We have checked that our result for iV = 2 is not changed if we take the values of a and ly taken in [Q. Also, the 
result does not depend sensibly on u*F^, although our value differs from their. We have traced back the difference 
between our result and |3^] to uB at the critical point: limiting ourselves to = 2, we have u*B* « 0.0391089 while 
p9| gives a value u*B* w 0.0363919. This difference is all is needed to explain the difference between our result and 



the result of Ref. l40f|, apart from a very small difference coming also from our use of Eq. (|104|) instead of (103). Since 



u* B* has been obtained in |39| using a five-loop Borel resummation, it is tempting to believe it is more accurate. For 
this reason, we have also determined numerically the five- loop strong-coupling limit of uB. We shall show a detailed 
numerical resummation in |^6|, showing here only the main steps. Starting from the five-loop expansion [ ^ , ^ 

uBiu) . + m±A^^ + A^(A^ + 2)(A. + 8)[-25+12C(3)]^, ^ ^ 
^ ^ 2 48 648 

(-319Ar2 _^ i3968iV + 64864 + 16(37V2 _ 382iv - 1700)C(3)96(4Ar2 ^ 39^ ^ 146)C(4) - 1024 (5iV + 22)C(5)) 5 

^ 41472 " ' 
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(153) 



and using the algorithm given by Eq. (O), the corresponding strong-couphng hmit is: 



N 





1 


2 


3 


4 


u*B* (5 loops) 





0.0209552 


0.0372717 


0.0502225 


0.0605918 


u*B* (Ref. Hi) 





0.020297 


0.0363919 


0.049312 





Our five-loop result is now much nearer to the Borel resummed values of | |39| than our three loop order of Section IV D . 
For this reason, we believe our five-loop result is near the infinite-loop limit extrapolation. More details will be given 
in j5^], which also contains the effect of variations of P'^, which is the second source, after u*B* , of error for Rc- 
Finally, our best values for the ratio Rc are collected in the next table: 



N 





1 


2 


3 


4 


Rc 





0.05803 


0.12428 


0.19402 


0.26285 


Rc (Ref. [|9)) 







0.123 


0.189 





To our knowledge no experimental value of this ratio is known for N = 2. The case TV = 3 is presented in Table 7.6 
of Ref. II^. For = 1, the value of the ratio has only slightly changed compared to the results based on three- loop 
a and v. This is due to the fact that, for = 1, a is positive and its effect on Rc is less sensitive. In the work [ p^ , 
the theoretical and experimental values of Rc are also given for = 1. The theoretical values seem to prefer a value 
around 0.057 while the experimental values are around 0.050. From Table 7.1 of |Q, we however see that values close 
to 0.06 might as well be obtained. 

Better experiments or other theoretical studies are needed in order to see if our predictions are correct or have to 
be ruled out. 

Finally, we conclude this section with the ratio of the susceptibilities for the Ising model. Using the critical exponent 
V to five loops (|l5l|), we obtain 



r+/r_ = 4.06419, 



(154) 



where we took the ratio (fx- / fx+T from (147). We might have slightly increased r+/r_ using the value 2.1447 



of Eq. (146). However, we would still be far from the value 4.77 of ||18|. The only possible quantity we may still 
vary in the ratio ( |105| ) is P1^. Our three- loop value is 0.745874, while the five-loop result given in using Borel 
resummation is 0.7568. Using this value in our formula for the ratio, we find 



r+/r_ = 4.27154. 



(155) 
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The ratio of the susceptibihties depends sensitively on P^. We postpone to |^6| the apphcation of variational pertur- 



bation theory up to five loops for the resunimation of and r+/r_. We do not however expect a resummed P; 



+ 



different from |39|. For this reason, the ratio (155) is probably the best we can obtain. A ratio of 4.77 obtained in 



and references therein seems to be ruled out from our analysis. 



V. CONCLUSION 



In this paper, we have shown that variational strong-coupling theory |^,^ can be applied not only to critical 
exponents, but also to various amplitude ratios. We have focused on two- and three-loop results were analytical 
results for the amplitude functions are known for all N both above and below T^. Our results are analytical 

expressions, except in the last section were we used more information to find A-^-/A^, Rc and r_|_/r_. The results are 
quite sensitive to the precise value of the critical exponents. In addition, a five-loop evaluation of the rcnormalization 
constant B* was necessary. The ratio Rc was so sensitive to it that a three-loop calculation was not sufficient. The 
same remark holds for P^, which affects mainly r+/r_. A numerical study of the known five-loop amplitudes will 



be done in , which will contain refined results compared to Section IV F 

One interesting observation of our work is that we can evaluate series which have caused problem in previous Borel 
resummations when the expansion coefficients in term of the renormalized coupling constant are not alternating to 
low orders. For these functions, strong-coupling theory turned out to work well. 

Having obtained analytical expressions in N, we have shown that the coefficient of the series in the bare coupling 
constant may vanish and change sign. At the two-loop level, this lead to diverging results near certain value of N. We 
have seen that the problem disappears at the three-loop level, because of the interplay of the different coefficients of 
the series. We could show precisely how it works because all our results were analytical and not restricted to integer 
values of N. 

When using variational perturbation theory, nothing is known on the nature of the optimal variational parameter, 
which can be a minimum, a maximum, or a turning point. The analysis performed here should help to identify the 
correct (numerical) solution at higher-loop order. See for example the amplitude f*^ , for which it is not yet clear 



which of the solutions (141) or (142) has to be chosen. 



48 



ACKNOWLEDGMENTS 

We thank B. Kastening, F. Nogueira and A. Pelster for valuable discussions during the completion of this work and 
M. Bachmann for help speeding up the calculations. 

APPENDIX A: EXPONENT ui FROM STRONG-COUPLING THEORY 

The power p/q of the leading power behavior ff^'^ of a function Wl whose perturbative expansion has been given 



in (20) can be obtained taking the logarithmic derivative, giving (21). A subtlety arises for functions going to a 
constant in the strong-coupling limit. For such functions, p vanishes and the corresponding /* in Eq. ( p7| ) vanishes. 
Care has to be taken: the limit /* ^ is different from imposing /* = 0. In the former case, we can identify q (or lu) 
by matching the series to achieve /* = 0. Working directly with a series which has f*~Q implies a leading behavior 



p' /q = —Lu/e. The algorithm (17) serves then to identify the coefficient Cq of the right-hand-side of Eq. ([181). As an 



example how to use the series, let us derive the relation |g,gUl 

-^-l = 4^?i^. (Al) 
e a log u B 

The left-hand-side is of the type (|l8|), and the algorithm can be applied. Formula ( |Al[ ) follows directly from (|T^). 
Alternative derivation starts from (|2l]): If p/q is vanishing, this means that its series has a leading exponent p' /q = 
— oj/e, which we derive in the following manner. Start from formula ( pl| ) with its exponent p' /q which we know from 
the general behavior ( p^ ) with /* = p/q = 0, i.e., 

y ^ _w ^ d\og{p/q) ^^2) 
q e d log u b 

where p/q is not yet taken at its asymptotic zero value, but is given as the right-hand-side of (|2l]). It then follows 

UJ , fW'l W'r\ , W'l p 



Taking the limit ub oo, the term p/q vanishes by hypothesis, and we end up once more with formula (Al). 

Although the algorithm (|l^ cannot be applied directly for the right-hand-side of (|l]) if p/q is vanishing exactly 
but only in the limit p/q — )■ 0, we can nevertheless use a trick to circumvent this problem: If the series for Wl has a 
vanishing leading power p/q, then Wl/ub has a power p' /q = —1. This allows to deduce 
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This shows that the right-hand-side of IjQVi) can be used to reach the Umit 0. Then, lu can be extracted either from (k1 



or from (Al). It is also clear from the expression (A3) that the right-hand-side has to be resummed blockwise: we 
have to use the intermediate result p/q — before tempting to resum. Using a full rcsummation of the right-hand-side 
of the latter equation would lead to badly resummed results (although the underlying e expansion would be the same): 



It is necessary to use p/q ~ in Eq. (A3), and not its analytical form which would have been mixed up with the 
power series of Wl"/W[. 



APPENDIX B: FREE ENERGY TO THREE LOOPS 

From the model Hamiltonian (^), the analytical calculation of the Gibbs free energy rB{m'^,UB, Mb) near the 
coexistence curve below Tc and for Mg = {'!>%) = above has been obtained at the two-loop order in [Q and at 
the three-loop order in [|9|, thus extending the = 1 calculation of Rajantic Q. We write directly the three- loop 
result: 



Tb = ^m'lMl + ubMI + Y,T.T.^-^^'^''''^bik{io,N)(24uBf-'{M^ 



3 6-1 1 



(24Mi 



In 



roL 



1 k 



{2AuBf 



, (Bl) 



6=1 1=0 k=Q 

where tol = m'^ -I- \2uBMg is the longitudinal bare mass, the transverse one tqt = + 4ub-^^s being included 
in the parameter Q — t^t / ^ol • The nonanalyticity in the coupling constant is seen in the last term. 

The functions Fuk can be found in [ p9| , p| . Since we need them later on in this Appendix, we shall write the nonzero 
components: 



^300 



1 -I- (Af - l)u;3/2 



^ 384^ 



-2{N - l)w^/^ + {N^ - l)w\ , 

^ _ 1 /l 



1 + 2wi/2 



ZOOTT 

- isik^ (^^ + 24^4 - - 'i^-'-'^' + 2Ar - 6 + 81n ^±f^ 



-6N -9 + 4:{N + 1) In ^ + 8 In 



2 + 2wV^ 



+ w{N -1) 



^ 23^ {3 + (iV - 1) [l + (A^ + 2)^V2] I ^ 
-F310 = ^f^^ ^9^' - 18 + 108 Li2 (^-i^ - (A^ - l)|4t()-i/2 + 4A^ + 2 - (A^ + 2) 



(B2) 
(B3) 
(B4) 
(B5) 



(B6) 
(B7) 
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- 12 Lis ( I ) - 321n2-6(ln3)^ + [lOiV + 32 - 16(2Ar + 3) ln2 + 481n3 - 8(iV + 1) Inw] 



+ -w(84A^- 100- 1281n2 



J 4321n- - 324Li2 

1658887r3 3 



(B8) 



432ci - 277r^ 



(iV- 1)( lQw~^/'^ + 



37V + 14 



7r2 + 18(ln3)2 + 36Li2 7 



16 



C2+4Li2-2-2Li2 I -- 



13 

6 In 3 - In 2 - y I In 2 



1 28 

— + 16w^/2 [7 - iV + (iV + 1) ln(16w) + 21n2 - 61n3] 
o 



+ Aw < 4c2 - 12Ar 



224 



6 61n3-ln2 



16 
15 



In 2 + 12 



2Li2-2 - Li2 - 



1 



(B9) 



where, in the coefficients i^aio and F320 below Tc, terms of order 0(w'^/^, ur'/'^ Inw) have been neglected (the coefficients 
are calculated in the vicinity of the coexistence curve where an expansion with respect to w is justified). The constants 
ci and C2 have been defined in the main text, see Eqs. (^6|) and (]97|). 
Inverting the equation of state 

d ^ 



gives, in the limit hs — > 0, and to this order, the square of the magnetization: 



(BIO) 



Mi 



ub 



10 - iV + A{N - 1) ln3 - 2{N + 2) In 



(24wb)2 



^.|(-2m'^)-V2 



19207r3 



27367V - 5904 - 6480ci + 240(7V - l)c2 - (TBN^ - 5N + 



1260 



(iV- l)Li2 - +9Li2 -- 



960(7V- 1) 



2 Lis (-2)- Li2 - 



630(iV- l)(ln3)2 -481n2 10(iV - l)ln2-60(iV- 1) ln3 + 1117V - 561 



240(12iV - 57) In 3 - 1440(iV + 2) In 



-2m'i 



(BU) 



(24ub)2 j 

The logarithmic terms in are nonanalyticities which can be removed using the length ^_ instead of m'^ < 0, 
see Up to the three- loop order, the relation between ^„ and ra'\ < is 

'1385 



2m'% = C 



N + 2 , iV + 2, 

1 + UB^- ^{uB^-f 



108 



7V + 2 
1087r3 



41n(24ui3^_) 



1 



3(4387V + 4349) + 576(iV + 8) Lia -- + 48(7V + 8)Tr^ + 8(437V + 182) In - 



(B12) 



Using (B12) in (Bll) one obtains an analytic function of ub, from which one extract the amplitude of Eqs. ( |85| ) 
and (psf), after proper normalization with the help of Z^. This has been done in |4nE3| and will not be repeated here. 
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The equations we have quoted here are mentioned because we shah need them below for obtaining the ampHtude 
functions f-^^ and f^_ which, to our knowledge, have not been determined analytically within this model. 

APPENDIX C: THREE-LOOP AMPLITUDE FUNCTION OF THE ISOTROPIC SUSCEPTIBILITY 

ABOVE Tc 



By definition, the amplitude of the susceptibility above Tc is obtained from the susceptibility at zero momentum 
f^^ = C+X^^B' where the inverse susceptibility is given by the two-point function F^' at zero momentum. The 
correlation length above the critical temperature ^+ is defined as in Refs. p9|,^j44|: 



e+ = X+,B{'l)dx+lBi<})/dq' 



(CI) 



Combining with the definition of f^^ , we have 



rB ^ ^^+]b 

Jx+ a„9 



(2) 



(C2) 



The derivative of r^-* with respect to is needed. This is in contrast to Refs. where only the combi- 

nation (pl|) was needed. For this reason, the intermediate result leading to Eq. (|C^) was not published . Being 



needed to determine the ratio Rc in (104) and the ratio of the susceptibilities (105), we derive it in the following 



The two-point function can be written as F^^* = rp + (/^ — ^^(q, ro,UB), where the self-energy has the expansion 
Ss(g, ro, wb) = X]m=i(~^s)'"S^'' (q, rg). The two-loop results have first been given in Appendix A of Ref. Q], with 
the result, up to order q^: 



r(2) =q^- + ro- 4(N + 2)ADUBrl^^ + 8{N + 2fAl-^ 



32 



iN + 2) 
(47r3) 



27r 27r q^ 
D-3 ~ 27 n) 



(C3) 



The pole at D = 3 can be eliminated by subtraction, leading to the masses to^ and m'^. This is however of no 
concern here since we are interested in taking the derivative with respect to q^: 



dV 



(2) 



N + 2u'j 



(C4) 



277r2 ro ' 

For the three- loop expansion, one must calculate the diagrams in Appendix B of Again, we concentrate on the 
derivative of the susceptibility at zero momentum, focusing on the diagrammatic Eq. (B5) of | |4C[ |. The corresponding 
vacuum diagrams have been given by Rajantie in |3^, see in particular its Eqs. (15) and (25) and, taking the 
appropriate derivative with respect to the mass, we obtain the contribution of the three-loop diagrams: 
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(3) 



g2=0 



(iV + 2)2 {N + 2){N + l 



547r3 



-8 + Stt^ + 32 In ^ + 36 Li2 ( 



-3/2 



(C5) 



This has to be combined with (C4) to yield the expansion 



dr 



(2) 



= 1 



N + 2u% 



277r2 



(iV + 2)2 (iV + 2)(A^ + l 



-8 + 3tt^ + 32 hi - + 36 Lia - 



1 



3/2 ■ 



(C6) 



27^3 547^3 

Since there is no hnear term ub, the amphtude of the susceptibihty above Tc to three loops requires only the 
one-loop order of the correlation- length ^_|_, i.e., the one-loop order of the susceptibility. To three loops, the following 
expression 



,2 ^_ 



1 



iV + 2 ^ A^ + 2, ^ ,2 
UB^+ H 2~('"s4+) 



l + 21n(24«Be+) 



N + 2 



3(3iV-f 22) - 144(7V 



Li2 



12(A^ -I- 8)tt^ - 2(437V + 182) In - 



(C7) 



is foimd in the literature, see [M. This is the analogue of (B12) above Tc- At the one-loop level, there is no distinction 



between r-g and m^, and we identify tq — ^^'^[1 + {N + 2)ub^+/t^]- Together with Eq. (|C6|), we arrive at 



/- 



dr 



(2) 



X+ 



= 1 



N + 2 2 (iV + 2)(iV + 8) 

-{ub£.+) - - 



277r2 



547r3 



-8 + 37r^ + 321n- 



36 Li, - 



(C8) 



This has to be compared with the numerical coefRcients a'm of Table 2 in ||5^] . 

Having calculated the bare amplitude function f^^ , we can now turn to the normalized one f~)^_^ . Since the latter is 
related to a two-point function, the normalization factor is equal to the wave function renormalization constant Z^: 
/^^ = Zc/if^^, with Zcf) being supplied by (|l6|). Using the relation ( [l3| ) between ub and ub — 'i*B'?+/(47r), we arrive 
at the normalized amplitude function of the susceptibility above Tc, expressed in terms of the reduced bare coupling 
constant ub- 



/x+ = 1 - :^(iV + 2)u% - ^{N + 2)iN + 1 



3 /I 

-113+ 1277^ + 1281n- + 144 Li2 — 

4 V 3 



27' ' 27 

The corresponding expansion in terms of the renormalized coupling constant gives 



(C9) 



8 

l--{N + 2)u'--{N + 2){N + 8) 



-21-1- 127r'' -I- 128 In - + 144 Lis 
4 ^ 



(CIO) 



where we used (|14|). Contrary to (|C9|) which is well-behaved regarding strong-coupling theory, Eq. (|C10D , which 
coincide with the numerical coefficients Cm-* of Table 4 of js^ , is problematic when considering the Borel resummation 
scheme: All its coefficients are negative. For this reason, we have not been able to reproduce the Borel resummation 



made in Ref. |£3|. We shall however make, in the main text, a comparison between the strong-coupling limit of (C9) 
and the resummation performed in |53]. 



53 



APPENDIX D: THREE-LOOP AMPLITUDE FUNCTION OF THE iV = 1-SUSCEPTIBILITY BELOW Tc 



In Ref. 13^, the amplitude function of the susceptibiHty below Tc for A'^ = 1 has been calculated numerically to 
five loops. We have quoted in (^Tj) the corresponding two-loop part. This amplitude function enters the ratio of the 



susceptibilities (lOf:). Since /^^ has been obtained to three loops in the previous section, it is also interesting to 



obtain analytically: The ratio ( |105| ) will thus be analytical. 

In Ref. 1^^, the free energy Tb has been given analytically up to three loops. We shall use this knowledge to 
determine f^- ■ We have recalled the relevant equations in the first part of this Appendix, which have to be evaluated 
for = 1 and u) ^ Q. The derivative of the free energy with respect to the magnetization leads to the equation of 



state (BIO) which can be inverted to obtain the magnetization . We have recalled its expression in (Bll). The 
equation of state can itself be derived with respect to the magnetization, defining the inverse susceptibility below 
Tc'- XZ^B ~ dhB /OMb- Only at this stage is the external field hB taken to be vanishing. The length ^_ [Q, which 
we recalled in ( p312| ), is then used to remove the nonanalyticity coming from logarithms of the coupling constant. 
Doing so, and using the magnetization given by the equation of state, we have been able to obtain the inverse bare 
susceptibility below Tc: xl s = ^-^/;f_i "^ith 



^ 2tt 367r^ 



19472- 11664ci + 37r^ + 10480 In - + 36 Lis ( 



647r3 



(Dl) 



Numerically, this expansion reads 



1 + 1.43239(uB^_) - 2.98616(wbC- 



11.2134(ubC- 



This result agrees perfectly with the numerical expansion given in the last column of Tabic 2 in 



(D2) 



Using the 



relation (14) between ub and ub, 'we obtain 

4244 , 



e_ = 1 + IMb - 



9 



19472 - 11664ci + Sn'' + 10480 In - + 36 Lia - 



1 



(D3) 



The renormalized version of ( p33| ) is found by multiplying it with from ([l^): 



/x_ = 1 + 18uB - 



4352 



19904 - 11664ci + Sn'^ + 10480 In - + 36 Lia - 



1 



(D4) 



This is the amplitude to be evaluated in the strong-coupling limit and entering the ratio ( 105 ) 



The bare amplitudes (C8) and (D3) might as well be chosen to enter the amplitude ratio (105) since the renor- 
maliazation constant Zj, drops out, being the same above and below Tc- We have however chosen to work with the 



renormalized quantities (|C9|) and (D4) 
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For completeness, we also state the expansion of fx_ in terms of the renormalized coupling constant. Using (Hj 



and (D4), we obtain 

= 1 + I8u 



1480 



1072- 11664ci + 371^^ + 104801n^ + SGLia I 



(D5) 



3 V 3, 

Taking the inverse of this equation, we recover the coefficients of the second column of Table 3 of Ref. . 

For an application to the evaluation of the amplitude ratio of the susceptibilities, we also give the perturbative 



expansion of the ratio fx-/ fx+ at iV = 1. Combining (C9) with (D4), we obtain 



fx- 



1 + ISUB 



1420. 



19184 - 11664ci + 997r^ + 9456 In - + 1188 Lia 

3 



(D6) 



APPENDIX E: AMPLITUDE RATIOS 



To get the different amplitude ratios of Section [VB, we make use of the relations 



X+ = Z^- — exp 



X- = Z^- — exp 

3X- 



('+: 



70 , , 



-y^ exp 



(El) 
(E2) 
(E3) 
(E4) 



which were derived in g| for Q, |38| for (|E|) and (^, and [|6| for (^. All the quantities have been defined in 
the main text, except for 

' r('±) du'\ 



u 



exp 



(E5) 



with u(l) = l± and the flow parameter chosen as l±fJ.£,± = 1, and with ^± = ^^\t 



Ij-l-v 



The amplitude ratio of the heat capacity (100) follows trivially from (E3), while the amplitude ratio for the 
susceptibilities (105) is a direct consequence of (El) and (|E2|). The only missing information is the ratio ^^J./^" , given 



explicitly in [R6| as 



e 



(E6) 



Because our derivation (104) of the universal combination Rc does not coincide with (103) derived by the authors 
of |40|, we reproduce below our calculation. We need the amplitude Am, related to (E4) by [|4j {(f)B) = Mb ~ AM\t\^- 
We deduce 
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70 



(E7) 



where we have specified that the right-hand-side is evaluated at the critical point. 



In the same way, the amplitude of the susceptibility is obtained from ( |El| ) using ||5J] x+ ~ r+|i| We deduce 

(E8) 



2± 



Taking the ratio /A\.j, we have directly 



A2 f* f* VS.- 

The dependence in |i| has disappeared, as it should, due to the identity — 2u + 2(5 — v{D ~2) — Q. 



Combining with the definition of in (E3), we get 



Rc 



V+A+ _ {h+f 



X (^-2) 



A 



(^+)2-.(g-2) 



f* f* 



(E9) 



(ElO) 



where we used (E6) to obtain the last equality. Using 6+ = 2uP+ and b = 3/2 — 2iyP^ pq], as well as A3 — l/(47r) 



from M), we arrive to the amplitude ratio Rc given in (104) 



APPENDIX F: DETERMINATION OF THE POLYNOMIAL P, TO THREE LOOPS 



In this section, we want to derive the analytical expression for the polynomial P+ up to three loops. It has been 
given numerically, and resummed, for = 1,2,3 up to five loops in ]5^ , so that our analytical result will have to 
match this reference. Above Tc, the relation between m'^ and the correlation length has been given in (C7) at the 
three-loop level. A polynomial Pf in powers of is defined through the relation 



pB =dm'l/di 



+ ' 



(Fl) 



leading to 

Pf = 1 
+ 



(UbU) n—iUBt.+ ) 



2lT 

N + 2 
1087r3 



-3{3N + 22) + 12{N + 8)tt'^ + 2{43N + 182) In - + 144(7V + 8) Lis ( -- 



(F2) 



4 ' ■ 'V 3, 

The numeric coefficient bm of Table 2 of | |53| coincide perfectly, up to three loops, with our analytical expression, 
which has the advantage of being valid for all N. Its renormalized counter part is defined by 



P+ = z-'p^, 



(F3) 
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where the renormahzation constant has been given to three loops in Eq. (p^). The corresponding power series 
in ub follows readily: 



^+ = 1 - 2(iV + 2)ub + 4(iV + 2)(2iV + 17)u 



3 /I 

-3(36iV2 + 8377V + 3920) + 2^(N + 8)'k^ + 4(43iV + 182) In - + 288(Ar + 8) Lis ( -- 



u\, (F4) 



where we have used the relation between ub and ub given in Eq. (13|), the scale /i being identified with the inverse 



of the correlation length: ub — Eq. (F4) is the polynomial whose strong-coupling expansion has to be 

calculated. The corresponding power series in the renormalized coupling constant u follows from Eq. (^): 



l-2(7V + 2)u + 4(7V + 2)u2 



+ ^(A^ + 2) 



-3(63/V + 572) + 24(iV + %)tx'^ + 4(437V + 182) In - + 288(iV + 8) Lis 



(F5) 



The reader can verify that the analytical result coincide, for TV = 1,2,3, with the numerical values in Table 4 of 
Ref. p3|. It differs only in the fifth decimal place of the cubic term cpz of this table. 
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